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Abstract 

The dynamics of cosmological models with isotropic matter sources (perfect fluids) is ex- 
tensively studied in the literature; in comparison, the dynamics of cosmological models with 
anisotropic matter sources is not. In this paper we consider spatially homogeneous locally 
rotationally symmetric solutions of the Einstein equations with a large class of anisotropic 
matter models including collisionless matter (Vlasov), elastic matter, and magnetic fields. 
The dynamics of models of Bianchi types I, II, and IX are completely described; the two most 
striking results are the following: (i) There exist matter models, compatible with the standard 
energy conditions, such that solutions of Bianchi type IX (closed cosmologies) need not nec- 
essarily recollapse; there is an open set of forever expanding solutions, (ii) Generic type IX 
solutions associated with a matter model like Vlasov matter exhibit oscillatory behavior to- 
ward the initial singularity. This behavior differs significantly from that of vacuum/perfect 
fluid cosmologies; hence "matter matters". Finally, we indicate that our methods can proba- 
bly be extended to treat a number of open problems, in particular, the dynamics of Bianchi 
type VIII and Kantowski-Sachs solutions. 



1 Introduction 

Spatially homogeneous (SH) cosmologies are one of the main cornerstones of General Relativity. 
On the one hand, they provide a rich supply of case studies for the effects that solutions of the 
Einstein equations of General Relativity can (and will) produce; in particular, the behavior of 
SH cosmologies gives the necessary input to model the behavior of our actual universe. On the 
other hand, there are convincing arguments that the dynamics of generic solutions of the Einstein 
equations close to spacelike singularities are built on the dynamics of SH cosmological models [3] . 

An SH cosmological model is a solution of the Einstein equations that is independent of the spatial 
variables in a suitable frame (see Sec. 2 for a precise definition). Therefore, for these models, the 
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Einstein equations reduce to an autonomous system of non-linear ordinary differential equations. 
As first recognized by Collins [10], provided the system possesses the appropriate regularity, the 
equations can be analyzed using the potent methods of dynamical systems theory. Since then, our 
knowledge of the (qualitative) dynamics of SH cosmological models has increased substantially, 
see [9, 17, 33] for reviews. 

Until recent years, the analysis of SH cosmologies has been largely restricted to isotropic matter 
sources like perfect fluids (including dust) and scalar fields. However, the results available for 
anisotropic matter sources, such as electromagnetic fields [20], collisionless (Vlasov) matter [15], 
elastic materials [4] , and viscous fluids [30] , suggest that the dynamics of SH cosmological models 
are in fact sensitive to the choice of matter model, not merely quantitatively, but also qualitatively. 
A first step towards a systematic analysis of the influence of the matter model on the (asymptotic) 
dynamics of solutions was taken in [5]: By analyzing Bianchi type I solutions (see Sec. 2 for the 
definition) of the Einstein-matter equations for a large class of (anisotropic) matter models that 
includes perfect fluids as a special case, it was shown that the qualitative dynamics of solutions 
(and the asymptotics towards the initial singularity, in particular) strongly depends on specific 
properties of the matter source. Differences between perfect fluid solutions and anisotropic matter 
solutions even exist when the matter model resembles a perfect fluid very closely [6] , which suggests 
that, at least within the family of matter models considered, the perfect fluid model is 'structurally 
unstable'. 

The purpose of this paper is to extend the analysis of [5] to the higher Bianchi types of class A. 
SH cosmologies of these types, especially the type IX solutions, are of particular interest. The 
important question of whether a cosmological model will expand forever or will recontract (and 
eventually recollapse) arises within this class of models. Moreover, the understanding of type IX 
solutions with their oscillatory behavior (Mixmaster behavior) toward the initial singularity is 
regarded as one of the keys towards a better understanding of singularities in General Relativity 
in general. 

In this work we restrict to SH cosmologies that are locally rotationally symmetric (LRS); in par- 
ticular, the isometry group of the spacetime is four-dimensional, see Appendix B. The main 
result of the present paper can be summarized as follows: The folklore statement "matter does 
not matter" has clear limitations, which can be specified rigorously. The qualitative dynamics of 
(SH LRS) cosmological models is in fact matter-dependent; cosmological models with anisotropic 
matter satisfying specific properties exhibit a completely different qualitative behavior than the 
corresponding perfect fluid solutions. The two most striking facts that we derive are the following. 

(i) There exist physically viable matter models such that the behavior of generic solutions toward 
the initial singularity is oscillatory; this is in stark contrast to the behavior of vacuum and 
perfect fluid solutions. These models arc LRS Bianchi type II and IX models. 

(ii) There exist physically viable matter models such that "closed- universe-recollapse" [2,21] does 
not hold; LRS Bianchi type IX solutions (closed cosmologies) need not necessarily recollapse; 
there exists an open set of initial data such that the associated solutions expand forever. 

By 'physically viable' matter models we mean matter models that satisfy the standard energy 
conditions, the weak energy condition, the dominant energy condition, and the strong energy 
condition in particular. Note that perfect fluid solutions of Bianchi type IX necessarily recollapse 
provided the strong energy condition holds. 

It is important to note that there exist important explicit examples of matter models that are 
directly covered by our analysis, e.g., elastic matter and collisionless (Vlasov) matter with mass- 
less particles, and matter models to which our techniques extend in a straightforward way, e.g., 
magnetic fields and Vlasov matter with massive particles; we refer to Sec. 12. 
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This paper is organized as follows. In Sec. 2 we define the class of spatially homogeneous space- 
times and the various Bianchi types; the Einstein equations for these models are given. In Sec. 3 
we specify the class of matter models we want to consider by mild and physically motivated as- 
sumptions. The family of matter models thus defined is sufficiently general to include important 
explicit examples and provides a non-trivial generalization of the perfect fluid matter model. (In 
Sec. 13 we elaborate on possible generalizations of our assumptions on the matter model.) In 
Sec. 4 we impose the condition of locally rotational symmetry (LRS) and discuss the Einstein 
equations and the matter variables for this case. In Sec. 5 we review the results on the dynamics 
of vacuum and perfect fluid models; the behavior of solutions with anisotropic matter we derive 
in the subsequent sections can be gauged against these results. In Sec. 6 we introduce a new 
set of dynamical variables, which recast the Einstein equations into an autonomous dynamical 
system — the reduced dynamical system — over a relatively compact state space; we remark that our 
variables are different from the standard Hubble normalized variables [33]. The following sections 
are devoted to the qualitative analysis of the reduced dynamical system for the different Bianchi 
types employing methods and concepts from the theory of dynamical systems. The lower Bianchi 
types, i.e., I and II, are studied in Sees. 7 and 8, respectively. The analysis of the Bianchi type IX 
case is divided into three sections: In Sec. 9 we slightly reformulate the equations for later pur- 
poses; Sec. 10 contains the analysis of the system, where we demonstrate the need for an additional 
variable transformation; finally, in Sec. 11 we collect the results and prove the main theorems on 
the dynamics of LRS Bianchi type IX solutions with anisotropic matter; in addition we put the 
theorems into a broader context. In Sec. 12 we show that our results directly apply to a specific 
matter model, namely collisionless (Vlasov) matter for massless particles; we thereby extend the 
results of [28]. Moreover, we show that our analysis extends straightforwardly to the magnetic field 
matter model. Finally, in Sec. 13 we discuss possible generalizations of our results and give a list 
of interesting open problems. The paper contains three appendixes. In Appendix A we present 
the exact solutions of the Einstein-matter equations that have been discovered in this work, whose 
role is that of 'attractors' of typical solutions. Appendix B contains a thorough discussion of local 
rotational symmetry. Finally, Appendix C reviews a few facts about dynamical system theory. 



2 Spatially Homogeneous spacetimes 

A spacetime (M, g) is spatially homogeneous if it admits an isometry group whose orbits are 
spacelike hypersurfaces that foliate M. The spatially homogeneous spacetimes divide into the 
Kantowski-Sachs models, see Appendix B.2, and the Bianchi models. 

A Bianchi model is defined to be a spatially homogeneous spacetime whose isometry group possesses 
a three-dimensional subgroup Q, dimC? = 3, that acts simply transitively on the spacelike orbits. 

Remark. There are three subcases. First, the Friedmann-Lemaitre(-Robertson- Walker) mod- 
els, see, e.g., [33, Chapter 2], which admit a six-dimensional isometry group. Second, the lo- 
cally rotationally symmetric (LRS) Bianchi models, see Appendix B.l, whose isometry group is 
four-dimensional; by definition, the group orbits coincide with the spacelike orbits of the three- 
dimensional subgroup Q. Finally, there are the Bianchi models that do not admit additional 
continuous symmetries beyond the three-parameter symmetry group Q. 

By definition, a Bianchi spacetime (M, 4 g) has the manifold structure M = I x Q, where / C M. is 
an interval. Let Gt, t € I, denote the Lie group Q t = {t} x Q, and let ^ = £i(t) be a left-invariant 
frame on Q t ; Latin indices run from 1 to 3. Expressed in terms of the coframe lo 1 = io l {t) that is 
dual to £j, the metric on Q t (the 'spatial metric') takes the form gij(t) uj 1 ^ . We temporarily adopt 
the Einstein summation convention, i.e., summation over repeated indices is understood. In the 
following we assume that the spacetime coordinate t is the standard 'cosmological time' ('Gaussian 
time'), i.e., the proper time along the geodesic congruence orthogonal to the group orbits. Then 
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the metric on the Bianchi spacetime (M, 4 g) takes the form 

4 g = -dt 2 +g %3 {t) wV» 



(1) 



The Einstein equations (in units such that 8ttG = c = 1) split into the constraint equations 

R - k ij k VJ + (g ij k %] ) 2 = 2T 00 , - V, (g hn k lm ) = -T 0j , (2a) 

which are called the Hamiltonian and the momentum constraint, respectively, and the evolution 
equations 

d tgij = -2kij , dtk\ = R i j + (g lm k lm ) ky - T) + (g lm T lm - Too) . (2b) 

Here, g^ = gij(t) is the (Ricmannian) spatial metric, g lJ its inverse; Rij is the associated Ricci 
tensor; R = g lm Ri m is the Ricci scalar; kij = hj(t) is the second fundamental form of the 
hypersurface of spatial homogeneity Qt', finally, Too, Th, and T^ are the energy density, momentum 
density, and stress tensor, which are derived from the energy-momentum tensor T^ v of the matter; 
Greek indices run from to 3. Above and in the following, (Latin) indices are raised and lowered 
with the spatial metric g^ and its inverse g lJ . 

We denote by C\j = C • (t) the structure constants of the Lie algebra associated with Q t . In this 
paper we consider Bianchi class A models, which are characterized by C j k = 0. (However, we also 
briefly discuss LRS Bianchi type III models which are of class B, see Appendix B.4.) Then the 
Ricci tensor Rij is given by 

Rij = -h&ti (C % + m m9 kn C m ni ) - \C m nk C p ql9jm g ip g kc >g ln . (2c) 

Let e^k denote the standard permutation symbol and define, as usual, n !J = J kl . Since 
Siji n kl = C k i3 + 28 k , i C^ m , it follows that, for class A models, 

C k ij =e i j l n kl , (3) 

where n kl is symmetric because £ m ij n ^ = C l ml = 0. Using (3), the momentum constraint in (2a) 
can be rewritten as 

e m n kl k\ = -T j . (4) 

The above equations are valid in a general left-invariant frame £j and its dual ui' 1 . The freedom 
of redefining such a frame by a t-dependent linear transformation can be used to simplify the 
equations, e.g., by introducing an orthonormal group- invariant frame (which is the approach used, 
e.g., in [33]). For our purposes it is preferable to exploit this freedom to make the frame time- 
independent. If the frame = £i (and w l (t) = u> 1 ) is time- independent, so are the structure 
constants C^- (i) = C* . and the matrix n y (t) = n 1 ^ ; we use the convention that hatted quantities 
are time- independent. Moreover, we are still allowed to use a time-independent (unimodular) linear 
transformation n tJ i— > U - k L J \n kl to diagonalize the constant symmetric matrix n y , i.e., without 
loss of generality, 

n lJ = diag(ni,n 2 ,n 3 ) . (5) 

After a further rescaling of the frame and up to a permutation or an overall change of sign, the 
('structure') constants hi, hi, h% are given as in Table 1. The values of these constants define the 
different Bianchi types of Bianchi class A. 

Remark. Using = C ij£k it is straightforward to see that the one- forms Co 1 satisfy 

dib 1 = —hi Co 1 A ill 3 , du) 2 = —hi a) 3 A a) 1 , dui 3 = —h% tu 1 A Co 2 . 

For each Bianchi type, the one-forms CS 1 admit different expressions in terms of local coordinates; 
we refer to [25]. 
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II 
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-1 


VIIo 
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VIII 


-1 


1 


1 


IX 


1 


1 


1 



Tabic 1: This table gives the classification of Bianchi class A models into the different Bianchi types. 
Each Bianchi type is defined by the values of the triple (ni,fi2, TI3); this is up to permutations and 
overall changes of sign of the constants (which correspond to discrete symmetries of the frame). 



We conclude this section by noting that the system of equations (2) implies a number of basic 
identities. Of particular relevance is the relation 

d t R = 2R ij k ij . (6) 

Since 1 

8 t R 2/,"/,',, • //TV/,, . 

to prove (6) we need to prove that g^dtRij = 0. Because we use a time-independent frame, 
the structure constants are time-independent; hence, differentiation of (2c), where we use the 
antisymmetry = —C^j, leads to 

gVdtRii = 9 11 C l kl C m nj {h m g kn - k kn g lm ) 

— I g u C T ^ lk C p ql ( — kj m gi p g kq g ln — ki p gj m g kq g ln + k kq gj m gi p g ln + k ln gj m gi p g kq } 

_ ril Am j. n ij n kn_s)l Am h.kn ij , Am Sip j, hq In _ rim rtP h.kq In 

— fci nj K im,g y o n jK g gi m + o „ fc o ql n pm g g o nfc o ql K g m pg 

= 0, 

which establishes the claim. 



3 Anisotropic matter models 

The system of equations (2) does not form a complete system in general unless some information 
on the energy-momentum tensor is given. For instance, one may specify the dependence of the 
energy-momentum tensor on some matter fields and add to the system (2) the evolution equations 
for these matter fields. These evolution equations are required to be compatible with the Bianchi 
identities V^T^ = 0; for certain matter models, such as perfect fluids, the Bianchi identities 
already contain the entire information on the dynamics of the matter fields. 

In this paper we do not explicitly specify matter fields or an energy-momentum tensor. Rather, 
the latter will be restricted by a set of mild and physically motivated assumptions which are known 
to hold for important examples of matter models. A particular consequence of our assumptions 
will be that the dynamics of the matter fields is contained in the Einstein equations (2) through 
the Bianchi identities. 
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3.1 Basic assumptions on the matter model 



The most common matter model in cosmology is a (non-tilted) perfect fluid. The energy-momentum 
tensor is 

r M v = p(dt® dt)^ + p[ 4 g + dt® dt] ^ , (7) 

where the energy density p and the pressure p are usually required to obey a linear equation of 
state, i.e., 

p = wp , w = const . 

The energy-momentum tensor (7) is isotropic in the sense that the eigenvalues of the spatial stress 
tensor T™. are all equal to the pressure p, since T l - = p8 l j- Moreover, the Bianchi identities 
V M T Ml/ = reduce to an ordinary differential equation for the energy density p, which can be 
solved to express p, and therefore the entire tensor in terms of the spatial metric: 

p = pon 1+w , p = constant , n = (det <?) -1 ' 2 . (8) 

In this paper we consider anisotropic matter sources that naturally generalize perfect fluids. The 
energy-momentum tensor is not specified explicitly; instead, we introduce a set of assumptions 
on the matter models that lead to a class of anisotropic energy-momentum tensors naturally 
encompassing (7). We note that there exist standard matter models that satisfy our assumptions; 
among these are collisionless matter (Vlasov matter) for massless particles and elastic materials 
(for a wide variety of constitutive equations); we refer to the subsequent remarks and to [5, 8]. 

Our first (and main) assumption is that, in analogy with the perfect fluid case, the energy- 
momentum tensor is given as a function of the spatial metric; this might be possible only under 
certain restrictions on the initial data. 

Assumption 1. There exist initial data for the matter such that the components of the energy- 
momentum tensor are represented by smooth (at least C 1 ) functions of the metric gij. We assume 
that p — p{gij) is positive (as long as gij is non- degenerate) . 

Remark. By Assumption 1, the components of the energy-momentum tensor are represented by 
functions of the spatial metric; note that the choice of frame may affect the form of these functions. 
In addition, these functions might include a number of external parameters or external functions 
(which describe the properties of the matter) and the initial data of the matter fields. Assumption 1 
implies that the matter does not add any dynamical degrees of freedom to the problem. We refer 
to [32] for a discussion of the Hamiltonian structure of the matter models satisfying Assumption 1. 

In the Bianchi type 1 case we have shown in [5] that the spatial components of an energy-momentum 
tensor satisfying Assumption 1 arc completely determined by the energy density 

Too =p = p(gij) 

through the formula 

We claim that eq. (9) holds for every Bianchi class A model. 

The key to prove eq. (9) is to use relation (6). Employing (6) in the derivative (w.r.t 
first equation of (2a) and using (2b), we easily obtain 

ftp = V'Crv + a>). (io) 

On the other hand, since p = p{gij) we may use the first equation in (2b) to get 

dtP = -2k i j g jl -^. (10') 



(9) 

time) of the 
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Equating the r.h. sides of (10) and (10') gives a polynomial of k % j whose coefficients depend only 
on gij. Therefore, the coefficients must be equal, which establishes eq. (9). 

The second basic property of the energy momentum tensor (7) that we want to require is the 
absence of heat flux and viscosity terms. 

Assumption 2. We assume that there exist initial data for the matter such that Toj = . 

Remark. In the Bianchi type I case this is not an assumption but a consequence of the field 
equations, cf. (4) with n u = 0. 



3.2 Diagonal models 

By the previous assumptions, the set of equations (2) forms a complete system of ordinary dif- 
ferential equations for the functions gij(t) and k % j(t). Initial data at to > are given by gij (to) 
and k l j(to) and the initial values of the matter fields (which determine the function p(gij) and 
thus the energy- momentum tensor). It follows from Assumption 2 and eq. (4) that the matrices 
n y and k % j commute; it is thus possible to make a change of frame i— > L\£k to simultaneously 
diagonalize the structure constants matrix n y and the second fundamental form fc*j(io)> this is 
simply achieved by choosing eigenvectors of k l j(to) as the frame; since these eigenvectors are — or 
can be chosen to be — orthogonal, we can assume that (to) is diagonal as well. The so constructed 
(time- independent, left-invariant) frame will be called the initially orthogonal frame. 

To ensure that (gij,k l j) remain diagonal for all times, the energy-momentum tensor must be 
diagonal in the initially orthogonal frame. 

Assumption 3. We assume that there exist initial data for the matter such that T % j is diagonal 
in the initially orthogonal frame. 

Assumption 3 ensures that the off-diagonal components of the metric satisfy a homogeneous ODE; 
uniqueness of solutions of the evolution equations then implies that (gij,k l j) remain diagonal for 
all times; hence, the initially orthogonal frame is in fact an orthogonal frame for all times. We 
refer to solutions of this type as diagonal models and restrict our analysis to these models. 

For diagonal models, the Ricci tensor i?y of the spatial metric g^, cf. (2c), is diagonal as well. To 
neatly express Rij in terms of g^ we define 



9a 9 JJ 9 



33 rihk 



•I. u >!>■■> V 9 

Here and henceforth we denote by (ijk) a cyclic permutation of (123) (i.e., Sij k = 1). Moreover, 
the Einstein summation convention is suspended, i.e., there is no summation over repeated indices, 
unless stated otherwise-. 



Based on (2c) we find that the components of the Ricci tensor are given by 

R\ = \ ni m i - ("i m j ~ hkink) 2 ■ (Ha) 

The spatial curvature scalar is 

R = fiimi (h 3 mj + n k m k - \riimi) - \ (n^nij - n k m k ) 2 . (lib) 

For diagonal models the Einstein equations (2) simplify. The Hamiltonian constraint, cf. (2a), 
becomes 

(trfc) 2 - V {k\) 2 +R = 2p, (12a) 
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where tr fc = J^i an d R is given by (lib). The first evolution equation, cf. (2b), reads 

dtgu — —2 ga fe I i , or, equivalently, dtirii — rrii ( tr k — 2k l ^j . (12b) 
The second evolution equation, cf. (2b), takes the form 

d t k\ = R\ + (tr k) k\ -pi+l^Pl-p) , (12c) 
where R 1 i is given by (11a) and where we have set 

T) = diag(pi,p 2 ,P3) ■ 

The quantities pi are called the principal pressures. (Note that when pi = p 2 = P3 ~ P we recover 
the orthogonal frame components of the energy-momentum tensor of a perfect fluid, see (7).) 

3.3 Properties of anisotropic matter models 

The principal pressures Pi give rise to the average (or isotropic) pressure 

p = |(pi +P2+Pz) ■ 
We make the following simplifying assumption: 

Assumption 4. We suppose that the isotropic pressure and the density are proportional, i.e., we 
assume p = wp, w = const, where 

we (-|,l). (13) 

Remark. The cases w = 1 and w = — ^ lead to different dynamics, which we refrain from discussing 
here. We remark, however, that our results for Bianchi type I solutions, see Sec. 7, are valid for 
w G (—1,1). In Sec. 13 we outline how to treat the more general case when w ^ const. 

According to Assumption 4, the energy density and the isotropic pressure behave like those of 
a perfect fluid with a linear equation of state satisfying the strong and the dominant energy 
condition. It is natural to focus the attention to anisotropic matter sources of this kind, because 
they generalize the class of perfect fluid models commonly used in cosmology, see [33, 17] and the 
references therein. 

Recall that the energy density of perfect fluid matter is described in terms of the spatial metric 
by (8). This relation possesses a natural generalization for anisotropic matter models satisfying 
Assumption 1-4. Let again n = (det g)^ 1 ^ 2 = \/g u g 22 g 33 and define 

Si :=(<? 11 +5 22 + 5 33 rV i . (14) 

Note that s± + S2 + S3 = 1, thus only two of the variables (si, S2, S3) are independent. The domain 
of these variables is the interior of the triangle 

T = {0 < Si < 1 : si + s 2 + s 3 = 1} . 

By Assumption 4 and (9) we obtain that p = pig 11 , g 22 , g 33 ) must satisfy the equation 

E^ = K 1 +->- ( 15 ) 

In the variables (n, Si, S2, S3) (which are related to (g 11 , g 22 , g 33 ) by a one-to-one transformation) 
eq. (15) reads nd n p = (1 + w)p. Integration yields 

p = n 1+w ^(s 1 ,s 2 ,s 3 ), (16) 
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where ip is a function on T; wc require ip to be smooth (at least C 1 ). 

This relation generalizes (8) (because it reduces to (8) when ip = const). We see that the energy 
density of an anisotropic matter model satisfying Assumption 1-4 depends on the spatial metric 
in a particular manner. Since ip is unspecified, (16) covers a large variety of matter models. 

Remark. Note that ip can (and in general will) depend on the initial data of the matter fields. 
This is true even in the simplest of cases, the case of a perfect fluid, where ip = po, cf. (8), where 
po = const. 

We define the rescaled principal pressures Wi to be 




evidently, wi + 1x12 + W3 = 3w. In general, the rescaled (anisotropic) pressures are not constant. 
However, as a consequence of Assumption 4, we observe independence of n. More specifically, by 
using (16) in relation (9), we find that 

1 \ 1 of dl °g^ 9logip\ , . 

Wi = Wi(si, s 2 , s 3 ) = w + 2 — Si } . (17) 

V O log Si d logs;/ 

Like the function ip, the rescaled pressures Wi arc defined on the domain T. For our purposes 
it is essential, however, that the quantities Wi be well-behaved in the limit (si,S2> s 3) - * 97, see 
Sec. 3.4. 

The rescaled pressures permit a straightforward representation of the energy conditions [22] . As 
stated in Assumption 1, our basic assumption is p > 0. Then the weak energy condition is expressed 
in terms of the rescaled principal pressures as > —1 and the dominant energy condition as 
\wi(si, S2, S3) I < 1, V(si, S2, S3) € T and Vi = 1, 2, 3. The strong energy condition is satisfied if the 
weak energy condition holds and if, in addition, w > —1/3. We are mainly interested in matter 
models that satisfy the energy conditions; however, we shall also discuss the more general case. 



3.4 Asymptotic properties of anisotropic matter models 

The matter models we want to consider are those that exhibit suitable asymptotic properties of 
the rescaled pressures Wi. 

Assumption 5. We assume that the functions Wi(si, S2, S3) representing the rescaled principal 
pressures possess a regular (sufficiently smooth) extension from 7 to 7. Furthermore, we assume 
that there exists a constant V— such that 

Wi(si, S2, S3) = V- when s, = , (18) 

for arbitrary i. 

Remark. The assumption implies that u),(si, S2, S3) = v+ := 3w — 2w_ when s, = 1, for arbitrary i. 

Remark. Some comments on Assumption 5 are in order. The assumption that the quantities 
Wi(si, S2, S3) take limits when (s 1; S2, S3) converges to a point on 87 corresponds to the assumption 
that the rescaled principal pressures take finite values under extreme conditions, i.e., when one or 
more components of the spatial metric gij go to zero or infinity (which occurs close to a singularity) . 
The second part of Assumption 5 means that the rescaled principal pressure in the direction i 
becomes independent of the values of gjj and gkk in the limit gu — > 00 provided that gjj, gkk 
remain bounded. (The complementary statement is implicit in the assumption. The fact that Wi 
is well-defined for Si = 1 means that the normalized principal pressure in a direction i converges 
to a limit as ga — > 0, when gjj and gkk remain bounded from below.) 
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lype 


fj V dl Llt^o 




D_ 


/3 < —2 


/3 = -2, io > | 


C_ 


/Je (-2,-1) 


W >0,/3>max{-2,-i±H} 


B_ 


= -l 


w > 


A_ 


/3e (-1,0) 


^>-i,/3>max{-l,-i±^} 


A 


/3 = 


to > -1/3 


A+ 


/?e (0,1) 


to > -1/3 


B+ 


= 1 


w > -1/3 


c+ 


/3e(i,2) 


violated 




/3 > 2 


violated 



Table 2: Anisotropy classification of matter sources according to the value of the anisotropy pa- 
rameter /?. The value j3 = is exceptional — it corresponds to a matter source which behaves 
like a perfect fluid 'at the singularity'. In the second column we give the range of j3 and w that 
guarantees the validity of the strong and dominant energy condition, see [5] for more details. 



Remark. Note that a linear equation of state connecting the principal pressures and the density is 
necessary for self-similar solutions to exist [24] ; Assumption 5 is formulated to allow for asymptotic 
self-similarity A natural generalization of Assumption 5 is discussed in Sec. 13. 

A quantity that will turn out to be ubiquitous in our analysis is the 1 anisotropy parameter' (3. As 
in [5] we define 



The anisotropy parameter /? is a measure of the deviation of the matter model from isotropy 
under extreme conditions (this is because w_ denotes the rescaled pressure Wi when g it — > oo, see 
Assumption 5). 

The anisotropy parameter (3 gives rise to a classification of the anisotropic matter models into 
different types, see Table 2 and Fig. 1. It is shown in [5] that, in the Bianchi type I case, this 
anisotropy classification is intimately connected with the qualitative asymptotic dynamics of the 
associated cosmological models. In this paper we will see that this connection carries over to 
higher Bianchi types to a large extent (although we shall need to slightly refine the anisotropy 
classification). 

We make a final simplifying assumption on the anisotropic matter model: 

Assumption 6. We assume that either there exists a unique isotropic state of the matter or the 
matter model is a perfect fluid. 



Assumption 6 can be formalized straightforwardly. Uniqueness of an isotropic state means that if: 
has only one critical point in the interior of T, which is either a maximum or a minimum. By (17) 
this entails the existence of a unique triple (si, 82,83) <G 7 such that wi = W2 = W3 = w at 
(si, S2j S3). This state is then a state of either maximal or minimal energy (where, physically, the 
more important case is that of a minimum of the energy). Although it is physically reasonable, 
Assumption 6 can easily be relaxed to admit the existence of more than one isotropic state; for 
simplicity we refrain from doing so. 

To conclude this section we remark that there exist important examples of matter models that 
satisfy our assumptions in the Bianchi type I case: Collisionless (Vlasov) matter for massless 
particles and elastic materials (for a wide variety of constitutive equations); we refer to [5] for 
details. These matter models fulfill our assumptions for the remaining Bianchi types as well, at 
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Figure 1: Anisotropy classification according to the value of the anisotropy parameter (3. (The 
special case A 1 * is discussed in Sec. 10; its range depends on the value of w.) 

least within the symmetry class that will be studied in this paper (locally rotationally symmetric 
models, see Sec. 4). In particular, elastic matter falls into one of the + types, depending on the 
constitutive equation of the material, whereas collisionless matter (for particles with zero mass) is 
of type B_|_ , sec Sec. 12. If the assumption of local rotational symmetry is removed, it is not known 
how to resolve the field (Vlasov) equation for collisionless matter, except in the Bianchi type I 
case (see [23]). This fact precludes to prove that collisionless matter satisfy our assumptions in all 
spatially homogeneous models (see the remark after Assumption 1). Finally there exist important 
matter models that do not satisfy our assumptions directly but to which our techniques can be 
easily extended. Examples of such models are Vlasov matter for massive particles, see [8], and 
magnetic fields, see Sec 12. 

4 LRS models with anisotropic matter 

In this paper we restrict our attention to locally rotationally symmetric (LRS) models. As discussed 
in detail in Appendix B, there exists a one-dimensional isotropy group in LRS models that defines 
a plane of rotational symmetry. By using an adapted frame one can choose two of the components 
of the diagonal metric to be equal (these are the components of the metric induced on the plane 
of rotational symmetry). W.l.o.g. we make a choice of LRS frame such that 

.922 = 533 • (20) 

To analyze LRS models with anisotropic matter we require that the matter model be compatible 
with the assumption of LRS symmetry; this assumption replaces Assumption 3. 

Assumption 3'. We assume that there exist initial data for the matter such that T l j is LRS, i.e., 
T % a is diagonal and T\ = T\ , in the LRS frame. 

The LRS models with anisotropic matter satisfying Assumption 3' are a subclass of the diagonal 
models. These models are generated by the prescription of LRS initial data, i.e., 522(^0) = 533(^0), 
k 2 2 (to) = k 3 3 (to), and Assumption 3'. 

Remark. Assumption 3' is satisfied if there exists initial data for the matter such that the function 
ip{s\, S2, S3) is symmetric under the exchange of the second and the third variable. 

The Ricci tensor of the spatial metric of a Bianchi class A LRS model is given by (112), i.e., 

-I 100 O A A TOO /\ 

R 1 = h 1 m 1 , R 2 = n\U2 m\mi — ^ n 1 . (21) 

The Einstein equations (12) simplify. For convenience we introduce variables that are commonly 
used in cosmology We define the Hubble scalar H and the shear variable a + to be 

H := -§ trk = -i (k\ + 2k\) , cr + = \ tr k - k\ = ± (k\ - k 2 2 ) . (22) 



11 



Then the Hamiltonian constraint (12a) becomes 

H 2 + | hifijininij = -|p + j^nfm 2 + a\ . (23a) 

where we have inserted R = R 1 1 + 2R 2 2 and (21). The evolution equations d t gu = —2guk l i , 
cf. (12b), turn into 

d t m 1 =mi(-H -4a+) , d t m 2 = m 2 (-H + 2a+) . (23b) 
Recall that g 11 = m 2 , g 22 = m\m 2 , cf. (112b). Finally, the evolution equations (12c) take the form 
d t k\ = i n 2 m 2 + (tr k)k\ - Pl + ± ( Pl + 2p 2 - p) , 

(23c) 

d t k 2 2 = i [Zhih-jmim-j - n 2 m 2 ] + (tr k)k 2 2 - p 2 + \ (pi + 2p 2 - p) , 

or, equivalcntly, using the constraint and expressing these equations in terms of H and cr + , 

d t H = -H 2 - 2al - i (1 + 3^)^ , 

6 (23c') 
d t a+ = —3Ha + + ^ (nlml — hih 2 mim 2 ) + (w 2 — w)p . 

The system of equations (23) represents the Einstein-matter equations for LRS models. Due to 
local rotational symmetry, the functions describing the properties of the anisotropic matter model 
under consideration, see Subsccs. 3.3 and 3.4, simplify significantly. 

Since 522 = ff33 implies S2 = S3, only one of the variables (si, s 2l S3) of (14) is independent. We set 

si = s 3 = s , si = 1 - 2s , (24) 

which implies that < s < i; note that s = g 22 /(g 11 + 2g 22 ) = (2 + m 2 /mi)^ 1 . With slight abuse 
of notation, we define 

V(s) :=V^(l-2s,s,s). (25) 
Employing (17) and (25) to compute the rescalcd anisotropic pressures we obtain 

Wl l-2s, s, s = w-2 l-2s , f yv ; , w 2 l-2s,s,s = w 3 l-2s,s,s = w+ l-2s — f-^ , 

c( log s a log s 

(26) 

where we have used that (dip/ds 2 )(l — 2s, s, s) = (dtp/ds^)(l — 2s,s,s), which follows from the 
symmetry of ip(si, s 2 , S3) in (s 2 ,s 3 ). 

Let us define the anisotropy function u(s) through 

u>2(l — 2s, s, s) = ^3(1 — 2s, s, s) = u(s) , wi(l — 2s, s, s) = 3ui — 2m(s) . (27) 

The anisotropy function u(s) encodes the properties of the rescaled anisotropic pressures; by (26) 
we have 

t \ , fi o N dl °gV'(s) 

u(s) = w + (1 - 2s) — . 

a logs 

It is this function that enters the Einstein-matter equations (23); its main properties will turn out 
to determine the dynamics of solutions of (23). 

Since < s < i, a priori, u(s) is defined on the interval (0,^). However, by Assumption 5 it 
admits a regular extension to the closed interval [0, ^] . By (27) we have 

s = : u(0) = V- , wi = w+ = 3u> — 2v~ , W2 = W3 = V- , 

s = ^ : u{ 2 ) = — ^ — = — 2 ' Wl = v ~ ' W2 = w 3 = — 2 — = — — t ■ 
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Remark. Let us comment on the physical interpretation of u(s) and v± . By (27), u(s) is the rescaled 
principal pressure in the directions tangential to the plane of local rotational symmetry; therefore, 
i>_ is the rescaled principal pressure tangential to the plane of local rotational symmetry — and 
v+ is the rescaled principal pressure in the orthogonal direction — under the extreme conditions 
represented by s — ► (i.e., .922 = .933 — * °o while gn remains bounded). Analogously, under the 
extreme conditions represented by s — ► \ (which corresponds to g\\ — > 00 and 1722 = 333 remaining 
bounded) the rescaled principal pressure tangential to the plane of local rotational symmetry is 
the mean of v± and the rescaled principal pressure in the orthogonal direction is u_ . 

The anisotropy parameter (3, see (19), takes several equivalent forms 

/3 = 2^^^2 W - M( ° ) ^4^ ) '". (28) 
1 — w 1 — W 1—10 

By Assumption 6, ip(s) has a unique critical point s G (0, ^); this is either a maximum or a 
minimum, i.e., ip'(s) = and tp"(s) ^ 0. The point s corresponds to a unique isotropic state: 
u(s) = w and thus w\ = u>2 = 103 = w for s = s. At s = s we obtain 

as s=s as s=s 

Therefore, we find that the anisotropic matter models we consider fall into two main classes: 

+ cases (/9 > 0) : u(0) <w < u{\) , u{s) = w , u'(s) > , (29+) 
— cases (j3 < 0) : u(0) > w > u(|) , u(s) = w , it'(s) < . (29-) 

Clearly, the + cases comprise the cases A_|_, B_|_, C_|_, and D_|_, sec Table 2 and Fig. 1; the — 
cases comprise the cases A_, B_, C_, D_. For j3 = 0, i.e., in the exceptional case Ao, we have 
u(0) = w = u(^); we distinguish two subcases: 

A 0+ : 13 = 0, u'(s) > 0; A _ : = , u'{s) < . 



5 Vacuum and perfect fluid dynamics 

The system of equations (23) constitutes the Einstein-matter equations for LRS models with the 
class of anisotropic matter sources under consideration. Setting the energy density (and the prin- 
cipal pressures) to zero, i.e., p — and pi = Vi, the equations reduce to the (SH LRS) Einstein 
vacuum equations. Setting pi = P2 = P3 = P = wp (or w\ = W2 = W3 = w) yields the (SH LRS) 
Einstein-perfect fluid equations. In this section we review the results on the asymptotic behavior 
of solutions of the system (23) for vacuum and for the perfect fluid matter models. For a proof of 
these results we refer to [33]. 

There are a few known explicit solutions of the Einstein vacuum and perfect fluid systems that 
will be relevant for our analysis. Let a, b denote positive constants. The following solutions are 
vacuum solutions, i.e., p = and pi = Vi, and of Bianchi type I. 

• The flat LRS Kasner solution ( Taub solution) T is represented by the spatial metric 

.9ii = a t 2 , g 22 = .933 = b . (30) 

• The non-flat LRS Kasner solution Q is given by 

gn=at~ 2 / 3 , g 22 = 533 = bt^ . (31) 
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In contrast, the following solutions are perfect fluid solutions (where the fluid obeys a linear 
equation of state), i.e., Pi = P2 = P3 = P = wp. 

• The isotropic perfect fluid solution (Friedmann-Lemaitre-Robertson- Walker solution) F is of 
Bianchi type I (i.e., hi = Vi) and reads 

4 

.911 = .922 = ff33 = at 3 d+™) , (32) 

where w £ (— 1, 1] is permitted. The Hamiltonian constraint implies p ~ gfij^jp t~ 2 , which 
is in accordance with (8). 

• The Collins-Stewart perfect fluid solution C is of Bianchi type II, i.e., h% = 1, fi2 = ^3 = 0; 
we have 

l — w 3+w rr , , — r 

gu=at^, g 2 2 = g 3 3 = btAi+™) with y/E = V( b > ( 33 ) 

where w £ (— i, 1] is permitted. The energy density is p — j^pfyz^" 2 - 

Remark. All solutions listed above are self-similar, see [33]. Moreover the are all singular at time 
t = 0. 

It is well known that for all Bianchi types except IX, vacuum solutions and perfect fluid solutions 
of (23) with initial data at t = to > become singular in finite time in the past (w.l.o.g. we 
assume the singularity to be at t = 0) and are defined for all times t > to in the future. In the 
Bianchi type IX case, solutions of (23) are singular both in the past and in the future. After these 
preliminary remarks, let us review the asymptotics of solutions of the system (23) in the vacuum 
and perfect fluid case. 

LRS Bianchi type I and VIIq. In the vacuum case a solution of (23) is either a Taub solution T, 
see (30), or a non-flat LRS Kasner solution Q, see (31). In the perfect fluid case every solution 
except the isotropic solution F, see (32), is asymptotic to a vacuum solution toward the singularity 
(i.e., as t — ► 0), while toward the future (i.e., as t — ► 00) every solution is asymptotic to the FLRW 
solution F and thus isotropizes. 

LRS Bianchi type II. In the vacuum case, every solution is asymptotic to T toward the singu- 
larity and to Q toward the future. In the perfect fluid case, every solution is asymptotic to the 
Collins-Stewart solution C, see (33), as t — > 00, while toward the singularity generic solutions are 
asymptotic to T; there is also a family of non-generic solutions that correspond to a set of zero 
measure of initial data, which behave like F toward the singularity and thus exhibit an isotropic 
singularity. 

LRS Bianchi type VIII. In the vacuum case, all solutions are asymptotic to T toward the singularity; 
while toward the future they are asymptotic to a 'plane wave' solution. (The plane wave solution is 
another (SH LRS) vacuum solution of the Einstein equations, see eq. (9.7) in [33]; it is of type III.) 
In the perfect fluid case, the asymptotics of generic solutions are the same. 

LRS Bianchi type IX. Toward both the initial and the final singularity, vacuum and generic perfect 
fluid solutions are asymptotic to T. There exists also non-generic perfect fluid solutions which 
behave like C or F toward both singularities. 

Remark. An important feature of the above results is that the asymptotics toward the singularity 
of generic Bianchi LRS perfect fluid solutions is the same as for vacuum solutions. This fact is 
referred to by saying that perfect fluid matter "does not matter" when the singularity of a Bianchi 
LRS cosmological model is approached. 

The main purpose of this paper is to extend the results on vacuum/perfect fluid solutions to the 
class of (anisotropic) matter models specified by the assumptions in Sec. 3. We shall see that matter 
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models of type Ao (which includes perfect fluids) have the same (generic) asymptotic behavior as 
perfect fluids. (Matter of type A "does not matter" when a singularity is approached.) However 
the other types of matter models may exhibit asymptotic behavior that differs considerably from 
that of perfect fluid (and vacuum) models (even for arbitrarily small (negative) values of the 
anisotropy parameter (3.) 



6 The reduced dynamical system 

In this section we formulate the evolution equations and the constraint in terms of normalized 
variables. The benefit: In these variables some equations decouple and thus the essential dynamics 
arc described by a system of equations that is of lower dimension (the reduced dynamical system). 
Moreover, the self-similar solutions of the system (23) (which govern the asymptotic behavior of 
general solutions) are represented by fixed points of the reduced dynamical system. The local 
and global stability properties of the fixed points of the reduced dynamical system are therefore 
intimately connected with the asymptotic behavior of the cosmological models under study. 

We define the 'dominant variable' D by 



D = niU2 3 mi " 12 . (34) 

The constraint guarantees that the square root is real, cf. (23a). We employ D to introduce 
normalized variables according to 

Hd = ^, S + = 1J< M 1 = ^>0, M 2 = ^>0. (35a) 
In addition to (35a) we define a normalized energy density f2 by 



^ = ^2 > , (35b) 



and we replace the cosmological time t by a rescaled time variable r via 

did , . 

Tr = Ddf (35C) 

Henceforth, a prime denotes differentiation w.r.t. r. Expressed in these variables the Einstein 
evolution equations (23b) and (23c') split into a decoupled equation for £), 



D' = -D^H D (l + q) + S + (l-H 2 D )j , (36) 
and a system of coupled equations for the normalized variables (35a), 

H' D = -(l-H 2 D )(q-H D S + ), (37a) 
E+ = -(2 - q)H D Z+ + (1 - H%)E% - \(-h\ M 2 + h x h 2 M X M 2 ) + 3fl(u(s) - w) , (37b) 

M[ = Ah (qH D - 4S+ + (1 - , (37c) 

M' 2 = M 2 (qH D + 2S+ + (1 - #f>)£+) , (37d) 

The evolution equation for s is 

s' = -6s(l - 2s)S+ ; (37e) 



15 



we are free to add this equation to the system (37a)-(37d). The quantity q ('deceleration param- 
eter') is given by 

q = 2T? + + i(l + 3w)n , (38) 

and the normalized energy density, f2, can be expressed in terms of the variables (35a) by the 
Hamiltonian constraint, cf. (23a), which reads 

E% + ±h\ M\ - \h x h 2 M x Mi + fi + (1 - H D ) = 1 . (39) 

It is important to emphasize that the variables in (37) are not independent; there exist two con- 
straints: The variable s, which appears as the argument of the function u(s), is determined by 

S =(2 + -) ; (40a) 
recall that < s < \. The second constraint is 

l-H D = (h 1 h 2 )^, (40b) 

which follows directly from (34). Using (40b) the Hamiltonian constraint (39) simplifies to 

S 2 + + ^M 1 2 + fi = 1. (39') 

Remark. A normalization that is more common than (35) is the Hubble-normalization. Hubble- 
normalized variables are defined as in (35a), where D is replaced by H. Accordingly, H D = 1 and 
the system of equations is given by (36) and (37b)-(37e) with (39), where H D is set to 1. 

Remark. It is sometimes useful to consider the evolution equation for il as an auxiliary equation, 
Of =ft(2H D q + 2T, + (l-H D )-(l + 3w)H D -6Y, + (u(s)-w)y (41) 

The system of equations (37) describes the dynamics of LRS Bianchi class A models. The equations 
for Kantowski-Sachs models and Bianchi type III LRS models derive from this system as well; for 
a brief discussion we refer to Sec. 9.2. 

The system (37) contains redundant information; we will see that the essential dynamics of every 
model can be described by a subset of equations of (37), which differs according to the Bianchi 
type and which we call the reduced dynamical system. 

In Bianchi types I and II, eq. (40b) enforces H D = 1 since h\fi2 = 0; therefore, the variables are 
Hubble-normalized. We restrict our attention to expanding type I and type II models, which are 
characterized by Ho = 1- (The case Ho = — 1 corresponds to models that are forever contracting. 
The dynamics of these models is obtained by replacing r with — r, i.e., by inverting the flow in 
the dynamics of the expanding models.) Bianchi type I is characterized by fi\ = and fi2 = 0. 
Accordingly, the dependence of the r.h.s. of (37b) on Mi and M2 disappears and the equations (37c) 
and (37d) decouple. Therefore, in the Bianchi type I case, the dynamical system (37) reduces to 
a two-dimensional system consisting of Eqs. (37b) and (37c). Bianchi type II is characterized by 
hi = 1 and n 2 = 0. Accordingly, the equation for M 2 decouples and we obtain a three-dimensional 
system consisting of Eqs. (37b), (37c), and (37e). Bianchi type VIo is incompatible with the 
assumption of local rotational symmetry (see Appendix B); LRS Bianchi type VILj reduces to 
LRS type I. Bianchi type VIII is characterized by h\ = 1 and h\fi2 = — 1. From (40b) and 
the positivity of Mi and M 2 it follows that Hd > 1. (We restrict our attention to expanding 
cosmological models and thus discard the case Ho < — 1 which corresponds to forever contracting 
models.) Bianchi type IX, on the other hand, is characterized by h\ = 1 and h\fi2 = +1, and (40b) 
implies that Ho G (— 1, 1)- Therefore, models of Bianchi type IX that are expanding initially need 
not expand forever; an eventual rccollapsc is to be expected; wc refer to Sec. 10. 
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Through the present approach, the analysis of the (asymptotic) dynamics of solutions of the 
Einstein-matter equations has been reduced to the problem of finding the a- and w-limit sets 
of orbits of the reduced dynamical system (for each Bianchi type) . Naturally, we are particularly 
interested in the asymptotic dynamics of typical orbits. 

We say that an orbit of a dynamical system is typical if its a-limit (resp. w-limit) set lies in the 
past (resp. future) attractor. (Solutions of the Einstein equations that are associated with typical 
orbits will also be called typical.) Atypical solutions are non-generic; the a/w-limit set lies outside 
of the past /future attractor. 

The subsequent sections are devoted to an in-depth analysis of the dynamics of LRS Bianchi class A 
models. We begin with a study of the dynamics of models of the lower Bianchi types. The results 
on the lower Bianchi types are results in their own right; in addition, this analysis is essential for 
the study of the dynamics of Bianchi types VIII and IX. 

7 Bianchi type I 

In this section we study the dynamics of LRS models of Bianchi type I. (Recall that LRS Bianchi 
type VILj is the same as LRS Bianchi type I, see Appendix B.) According to Table 1, Bianchi 
type I corresponds to setting hi = and t%2 — fi3 = in (34) and in the reduced dynamical 
system (37). We see that Hd = 1, which reduces the normalization (35) to the standard Hubblc- 
normalization, and the equations for M\ and Mi decouple from (37). Therefore, the essential 
dynamics are described by the equations for the pair (£+, s), which are given by 

£'+ = -3(1 - Z 2 + ) - w)S+ - Ms) - w)] , (42a) 
s' = -6s(l - 2s)£+ . (42b) 

The Hamiltonian constraint (39') becomes 

z'l + n = i . (43) 

Therefore, the reduced dynamical system (42) is defined on the state space 

A' I = {(E + , S )e(-l,l)x(0,i)}, (44) 

which is a rectangle. 

The system (42) admits a regular extension to the boundary of the state space. We distinguish 
the vacuum boundary, which is given by = (i.e., p = 0), or, more specifically, by Vi = 7~U Q, 
where 

T = {£+ = -1} x (0,i) and Q = {£+ = 1} x (0, §) , 
and the matter boundary, which is the boundary subset satisfying S! =/= 0; it consists of 

I b = (-1,1) x {s = 0} and Zj, = (-l,l)x {«=§}. (45) 

The number of fixed points on the boundary of X\ varies between four and six, depending on the 
anisotropy parameter /3, see (28). 

T b /T tt There exist two 'Taub points': T b : (£+, s) = (-1, 0) and T s : (£+, s) = (-1, ±). 
Qb/Qu There exist two 'non-flat LRS points': : (S + ,s) = (1,0) and Qj : (£+,s) = (1, i). 
R^ The fixed point R^ : (£+, s) = (— /3,0) exists when (3 £ (—1, 1). 
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Rj The fixed point Rj : (£+, s) = (f , ±) exists when /3 € (-2, 2). 

The T-points are associated with the flat Taub solution (30), whereas the Q-points are associated 
with the non-flat LRS solution (31). The fixed points R,, Rt,, when they exist, are associated with 
non- vacuum solutions (since 17 > at these points), see Appendix A. 

Remark. We adhere to the convention that fixed points and invariant subsets with s = are 
denoted by a subscript b , while the subscript , indicates that s = i. 

Consider, first, the vacuum boundary Vi. The orbit T with S + = —1 connects the fixed point 
T\, with the point Tj. Like the fixed points themselves, the orbit T is a representation of a Taub 
solution (30) in terms of the normalized variables. Likewise, the orbit Q with S + = +1 connects 
the fixed point Q, with the point Q|, and represents a non-flat LRS solution (31). 

Second, consider the boundary components Jj, and Zjj. Using (28) we have 

l b : £; |s=o = -§(l-5^)(l-™)(£ + +/3) ! (46 b ) 

Z$ ■■ S+ \s=i/2 = -K 1 - - - f) ; (46,) 

in particular, R^ and Rj, when present, attract solutions on T\, and Zjj, cf. Fig. 2. To determine 
the role of the points R|, and R, for the flow on X\, i.e., for (42), we use that 

[f logs] |Rb = 6/3, [£bg(l-2*)] |R| =3/?. 

We infer that, in the — cases, i.e., if f3 < 0, these fixed points attract orbits from the interior of 
X\ and thus represent sinks in X\, while they are saddles in the + cases, i.e, if j3 > 0. (In the 
cases Ao± the fixed points are non-hyperbolic sinks and saddles, respectively, which follows from 
a center manifold analysis.) 

In the case D_|_ , the boundary of X\ forms a heteroclinic cycle connecting the fixed points at the 
vertexes of the rectangle, see Fig. 2. We schematically represent this heteroclinic cycle as 



i 


1 


T b «— 





We shall sec that heteroclinic cycles of the form (47) appear on the boundary of the state space in 
all Bianchi types. 

By Assumption 6 there is one and only one fixed point in the interior of X\. 

F The fixed point F is given by (S+, s) = (0, s), where s £ (0, h) is the unique solution 
of u(s) — w (<=>■ ip'(s) = 0). 

The fixed point F is associated with the isotropic perfect fluid solution (32). The eigenvalues of 
the linearization of the dynamical system (42) at F are given by 

A± = -| (l - w ± - w) 2 - 325(1 - 2s>'(s)) . 

Since u'{s) ^ tp"(s) 7^ 0) by (29), the fixed point F is hyperbolic. From (29) it follows that 
F is a sink for (3 > (i.e., in the + cases) and a saddle for (3 < (i.e., in the — cases). Moreover, 
in the former case, if u'(s) is large, A± is complex and the approach of solutions to F is oscillatory. 
(The cases Aq+ and Aq— behave like the + and — cases, respectively.) 
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Type 


Past Attractor 


Future Attractor 


D_ 


Tb, Qtt 


Tj, Qb 


C_,B_ 


Tb, Qfl 


Rj, Qb 


A_, A - 


T b, Qtt 


Rjj, Rj, 


A+, A 0+ 


T b, Qtt 


F 


B+C+ 


Q« 


F 


D+ 


Hcteroclinic cycle (47) 


F 



Tabic 3: The past and the future attractor for Bianchi type I orbits encode the asymptotic behavior 
of typical solutions. The possible a/w-limits of non-generic orbits can be read off from Fig. 2. 



Theorem 1. The phase portraits of the dynamical system (42) for the various anisotropy cases 
are correctly depicted in Fig. 2. The past and future attractors in the various cases are listed in 
Table 3. 

Proof. Consider the function 

Z 1 = (l-S2)- 1 V.(s). (48) 
A straightforward calculation shows that 

Z[ = -3(1 - w)Vl Z x , Z'(' |s+=0 = -54(1 - w)(u(s) - w) 2 ^(s) . 

It follows that the function Z\ is strictly monotonically decreasing along orbits in A"i\{F}. Applying 
the monotonicity principle, sec Appendix C, we infer that the a- and aj-limit sets of every orbit 
of the dynamical system (42) must be located on dX\ U {F}. Using the local analysis of the fixed 
point F and the straightforward dynamics on dX\ yields the claim. □ 



Interpretation of Theorem 1. In the + cases, the orbits of the dynamical system converge to the 
fixed point F. This means that each solution of the Einstein equations isotropizes toward the 
future, i.e., the metric is asymptotic to the isotropic perfect fluid solution (32) as t — -> oo. This 
is also the future asymptotic behavior of Bianchi type I cosmological models with perfect fluid 
matter, see Sec. 5. (Note, however, that isotropization rates may differ considerably [6].) In the — 
cases, in general, cosmological models do not isotropize toward the future (which is irrespective of 
whether the energy conditions are satisfied or not). In particular, in the cases B_, C_, D_, there 
exist typical solutions of the Einstein equations, see Sec. 6, that are future asymptotic to vacuum 
solutions. Other typical solutions approach the models corresponding to R|, and Rj; we refer to 
Appendix A where we explicitly give the metric associated with these models. 

Furthermore, the theorem shows that the past asymptotic behavior of models is intimately con- 
nected with the behavior of the vacuum solutions (30) and (31), which is analogous to the case of 
perfect fluid cosmologies. However, there is an interesting exception: In the case of a high degree 
of energy condition violation, i.e., case D_|_, the approach to the singularity is oscillatory. The 
solutions oscillate between the Taub family (30) and the non-flat LRS family (31) of solutions. 
Since Q, -/-> for these solutions, this is an example of asymptotic behavior where "matter mat- 
ters" , which means that the asymptotics are dissimilar from the asymptotics of vacuum models. 
For more details and for a generalization of the result of this section to the non-LRS case we refer 
to [5]. 

Remark. The LRS subcase of Bianchi type VIIo coincides with LRS Bianchi type I. Setting hi = 
in (37) yields the type I system (42) irrespective of the value of hi = hs. 
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Figure 2: Phase portraits of orbits representing LRS Bianchi type I models; the state space is 
X\ and the dynamical system is (42). Continuous lines represent typical orbits — each of these 
orbits is a member of a one-parameter set of orbits with the same asymptotic behavior. Dashed 
lines represent orbits whose future and/or past asymptotics are non-generic. Certain features in 
these diagrams are included for future purposes: First, the fixed points are color-coded. When 
the Bianchi type I state space X\ appears as the boundary subset of the type II, there exists an 
orthogonal direction. A black fixed point is an attractor, a white fixed point is a repellor in the 
orthogonal direction. Second, the cases A_ and A_|_ are split into subcases, which become relevant 
for Bianchi type II; note that the difference is merely in the color-coding (which is irrelevant for 
type I itself). 



8 Bianchi type II 

According to Table 1, Bianchi type II corresponds to setting hi = 1 and hi = h^ = in (34) and in 
the reduced dynamical system (37). We infer that Hp = 1, which reduces the normalization (35) 
to the standard Hubble-normalization, and the equation for M2 decouples from (37). Therefore, 
the essential dynamics are described by the equations for the variables (£+, Mi, s), which arc given 
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Figure 3: The state space X\\. The base of the 'tent' is the type I state space X\\ the 'roof Vn is 
the vacuum boundary; the sides are denoted by S\, and Sj, respectively. Recall that the subscript b 
corresponds to s taking the value s = 0, while jj corresponds to s = 1/2. 

by 

£' + = ±Mf (2 - £+) - [|(1 - - (u(s) - w)] , (49a) 

s' = -6s(l - 2.s)S+ , (49b) 

M( = Mi [2 (1 - 2E+) - |M| - |(1 - w)fi\ . (49c) 

The Hamiltonian constraint (39') becomes 

Y? + + i M x 2 + = 1 . (50) 

Therefore, the dynamical system (42) is defined on the state space 

#11 = {(£+,«, Mi) | -K £+ < 1, < s < \, 0< Mi < ^12(1 -£2) } ) ( 51 ) 
which resembles a tent; see Fig. 3. 

The dynamical system (49) admits a regular extension to the boundary of X\\, which consists of 
four components, 

dX n = XzUSt, US$ U Vn , (52) 

which are defined as 

Xi = {Mi = ; -1 < £+ < 1, 0<s<|} 

S b = {s = ; -1< £+ < 1, < Mj < ^12(1 -£2_) } , 

S t = {s = \ ; -1< £+ < 1, < Mj < ^12(1 -E^.) } , 

Vn = {O = ; -1< £+ < 1, < s < |, Mi = ^12(1 - E^_) } . 

The symbol X\ is used with slight abuse of notation; the type I state space X\ is embedded into 
X\\ as the boundary subset {Mi = 0}, see the discussion below. The component Vn is the vacuum 
boundary; the components Si, and Sj we call the 'sides'; see Fig. 3. Let us study the dynamics of 
the flow induced on each of the above boundary subsets. 
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The Bianchi type I boundary X \. The dynamical system induced by (49) on the invariant 
subset Mi = is identical to the type I system (42); its dynamics are described by Theorem 1, 
i.e., Fig. 2. We thus refer to X\ as the Bianchi type I boundary of X\\. Let us replicate the list of 
fixed points: 

T b /T tt The 'Taub points': T b : (£ + ,s,Mi) = (-1,0,0) and T tt : (£+,s,Mi) = (-1,±,0). 
Qb/Qtl The 'non-flat LRS points': Q b : (£+,s,Mi) = (1,0,0) and Qj : (E+,s,Mi) = (1, |,0). 

R b For p £ (-1, 1) we have R b : (£+, s, Mi) = (-/?, 0, 0). 

Rj For /3 G (-2, 2) we have Rj : (£+, s, Mi) = (f , ±, 0). 
F The Friedmann point F is given by (£+, s, Mi) = (0, s, 0), where s solves u(s) = w. 

It is straightforward to study the role of the fixed points w.r.t. the direction orthogonal to X\\ 
from (49) we have 

[^r(logMi)] {Xj = 2(1 - 2E+) - |(1 - w)(l - E* ) = i(i + 3w) + E+(f (1 - W )E + - 4) . 

As an example, consider the point F: Evaluated at F, this derivative is positive, so that there exists 
(at least) one orbit that emanates from F into the interior of X\\. In other words, F is a repellor 
in the orthogonal direction (i.e., the Mi direction); this behavior is color-coded into Fig. 2 — the 
fixed point F is represented by a white circle. (Points that act as attractors in the orthogonal 
direction are depicted in black in Fig. 2.) 

The vacuum boundary V n- The dynamical system induced on Vn is obtained by setting O = 
in (49); since O = entails p = 0, we refer to Vn as the vacuum boundary. The essential dynamics 
are described by the system 

E' + = 2(1 - E^)(2 - £+) , s' = -6s(l - 2s)£+ , 

since O = implies Mj 2 = 12(1 — E+). The closure of Vn is the rectangle (E + , s) e [—1, 1] x [0, i]; 
the vertices are the only fixed points: T[,, Th, Qb, Qj- The boundary of the rectangle is characterized 
by a simple flow, see Fig. 4. 

The orbits T[, — > Tj and Qi, <— Qj coincide with the type I vacuum orbits and thus represent the 
Taub solution (30) and the non-flat LRS solution (31), respectively Analogously, every orbit in 
the interior of Vn represents a type II vacuum solution. Since E' + > in Vn, the analysis of the 
dynamics is straightforward: The fixed points Tt, and Qi, are the a-limit set and the aj-limit set, 
respectively, of every orbit, see Fig. 4. 

The side S%. The dynamical system that is induced by (49) on <Sy is given by 

EV = ±Mf(2 - E+) - §0(1 - «;)(E + - |) , (53a) 
M[ = Ah [|(1 - w)E^ - 4E + + 1(1 + 3w)(l - ^Mj 2 )] , (53b) 

where = 1 — E^ — j^Mf. The state space is (the interior of) a closed semi-circle. The two 
vertices are the fixed points Tj and Qj. 

On the vacuum part O = of the boundary of <Sjj (which is the semi-circle) we observe 

S V iffco = 2 (! - S2 +X 2 - S +) ' ( 54 ) 

which is positive since |E+| < 1. Therefore, the semi-circle corresponds to an orbit Tj — > Qj; see 
Fig. 5. (The same orbit is depicted in Fig. 4 and the color-coding in Fig. 2 is consistent as well: 
In Fig. 2, Tj is white and Qj is black.) 
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Figure 4: The flow induced on Vn. 



The equation induced on the subset Mi = of the boundary of 5jj (which is the base of the 
semi-circle) is 

^ |Ml=0 = -§(i-£ 2 + )(i- W )(£ + -f), 

which coincides with (46|j) ; this is because the base of the semi-circle coincides with the Ij subset of 
dX\. Accordingly, there is a fixed point, Rjj, if \(3\ < 2 (i.e., in the cases A, B, C), which acts as an 
attractor on this boundary component; in contrast, we obtain an orbit Tj <— Qjj for (3 < —2 (D_) 
and an orbit Tj — > Qj for (3 > 2 (D_|_); see Figs. 2 and 5. 

To investigate which role the point Rjj plays in the context of the flow on <Sj we consider 

[£ logMx],^ = §/? 2 (l-™)-2/?+±(l + 3 W ). (55) 

Let 



8 - 2V9^ 2 - 6w + 13 8- 2^/(3(1 - w)] 2 + 4(1 + 3w) ^ , n ^ 

m w )~ ^71 \ = \ e (0,1). 56) 

3(1 — w) 3(1 — w) 

The r.h.s. of (55) is positive for —2 < (3 < [3$(w); hence, in this case, there exists one orbit that 
emanates from the fixed point Rjj into the interior of S$, i.e., Rj is a saddle in Sj. If, on the other 
hand, f3$(w) < (3 < 2, then the r.h.s. of (55) is negative, which makes Rj a sink in <Sj. (In the 
exceptional case (3 = f3$(w), the fixed point Rj is a non-hyperbolic sink, as will be proved below 
by using global methods.) Note that the distinction (3 ^ (3$(w) leads to two subcases of the case 
A_|_, since (3^{w) € (0, 1), cf. Fig. 5 and the color-coding in Fig. 2. 

Let us turn to the interior of St. For (3 < (3§{w) (i.e., in the — cases and the respective subcase of 
A_|_ ) there exists a fixed point in the interior of iSj , which we call Cj : 

C B The coordinates of C, are E + = j^^, = ^-^^Z^ 1 ^ ■ 

This fixed point represents a solution that generalizes the Collins-Stewart solution (33); it reduces 
to (33) in the case (3 = 0. The exact solution that corresponds to the fixed point Cj is calculated 
in Appendix A. The condition (3 < (3§{w) assures that (Mi) |c > 0; the point Cjj converges to Rj 
when f3 f3$(w). The point Cj is a sink in 5j; this will be proved below by using global methods. 
Finally, to investigate which role Cj plays in the context of the system (49) on the entire type II 
state space X\\, we compute 

[£log(l-2*)] |C8 =6S +|C9 >0, 
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which implies that Cj repels orbits from the interior of X\\\ in other words, Cj is a saddle in X\\. 
This fact is depicted in Fig. 5 by representing Cj by a white circle. 

To perform a global dynamical systems analysis of the flow on the invariant set Sj, we apply the 
monotonicity principle. 

Consider first the case D_|_. Since > 2 in this case, we have > in the interior of 5j. 
Therefore, the a- and w-limit set of every interior orbit must lie on the boundary dS$. From the 
simple structure of the flow on dS$, we deduce that the a- and w-limit set of every orbit in the 
interior of Sj is the fixed point Tj and Qj, respectively; see Fig. 5. 

In the remaining cases, where < 2, we define 

_ 2(l+3m) r _ 3(l-tu)(c-g/2)c 

C — 16-3^(1-10) ' ~ (l+3w)(l-c 2 ) ' 

It is elementary to verify that c £ (0, 1) and S < 1 for all < 2; moreover, 

8 < for £ (J3$(w), 2) , 5 = for = % (w) , and 6 > for < fy{w) . 
Consider the function 

2 ~ (l-cE + )2 ■ 

(For = the function Z 2 reduces to the monotone function for perfect fluid models given in [33, 
eq. (6.47)].) A straightforward calculation shows that 

Z'= X Z 2 , where x = f^^m (s+ _ c) 2 > Q (V/3 < 2) 

and % = if and only if E+ = c. If /3 < 0$(w), the fixed point Cj is located in the interior of <Sjj. 
In this case the function Z2 is monotonically increasing on <Sjj\{C|j}. The monotonicity principle 
then implies that the fixed point Cj is the w-limit set of every orbit in the interior of iSjj , whereas 
the a- limit set is contained on the boundary dS$. The generic a-limit set is the point Tj in the 
cases A + (where < 0^(w)), A_, B_, C_; see Fig. 5. In the case D_, the entire boundary dS$, 
which forms a heteroclinic cycle, is the a-limit set. Finally, for £ [0«(w), 2), which comprises 
the respective subcase of A+ and the cases B+, C_|_, the function Z2 is monotonically increasing 
in the whole interior and thus the a- and the w-limit sets of interior orbits are both contained on 
the boundary of S$. The a-limit set is Tjj, the w-limit set is Rj. These results are summarized in 
Fig. 5. 

The side S\,. The dynamical system induced on S], is given by 

S' + = \Ml{2 - S+) - §0(1 - «;)(£+ + 0) , (53a') 

M[ = Mi [f (1 - w)Y? + - 4S+ + |(1 + 3u;)(l - T^AIi)] , (53V) 

where CI = 1 — — j^Mf. The state space S\, is an identical copy of the state space Sj, where the 
dynamical system (57) on S\, is obtained from the system (53) on S$ upon the formal substitution 
1 — ► —2/3. Accordingly, we simply adapt the results of the analysis of <Sj by replacing 1— > —20. 

The flow on the vacuum part fl = of dS[, (which is the semi-circle) is determined by an equation 
identical to (54); hence, the semi-circle is an orbit T|, — > Q|>; see Fig. 6. The part Mi = of dS\, 
(which is the base of the semi-circle) coincides with the I\, subset of dX\] the equation is (46^). 
Accordingly there is the fixed point Rj, if \0\ < 1 (i.e., in the A cases), which acts as an attractor 
on this boundary component; if < — 1 (B_, C_, D_), we obtain an orbit T|, — > Qi,; if > 1 
(B_|_, C_|_, D_|_), we obtain an orbit Tj, <— Q;,; see Figs. 2 and 6. 

To investigate which role the point R[, plays in the context of the flow on Si, we consider 

[£ log Mi] |R = f/? 2 (l - w) +40 + |(1 + 3w) . (55') 
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(a) D_ 



(b) A _, A_, B_, C_ 




Figure 5: Phase portraits of the flows on <Sj in dependence on the anisotropy parameter [3. Con- 
tinuous lines represent typical orbits — each of these orbits is a member of a one-parameter set of 
orbits (i.e., a two-parameter set of solutions). Dashed lines represent orbits whose future and/or 
past asymptotic behavior is non-generic. The fixed points are color-coded. When the set 6>j i s 
viewed as the boundary subset of the type II state space Xn, there exists an orthogonal direction. 
A black fixed point is an attractor, a white fixed point is a repellor in the orthogonal direction, 
cf. Fig. 2. Analogously, the set Sj appears as a boundary subset of the type VIII and the type IX 
state spaces; the color-coding then refers to the orthogonal direction of when viewed as an 
embedded set in these state spaces. 



Let 

AH '= 3(1^) = 3(1^) £ ( "- 0) • (56 } 

The r.h.s. of (55') is positive for (3\,{w) < (3 < 1; hence, in this case, there exists one orbit that 
emanates from the fixed point Pq, into the interior of iSb, i.e., Rb is a saddle in S\,. If, on the other 
hand, — 1 < (3 < (3t,(w), then the r.h.s. of (55') is negative, which makes Rb a sink in S^. (The 
fixed point Rb attracts interior orbits in the case (3 = (3\,{w) as well, as will be proved below by 
using global methods.) Since (3\,{w) G (—^,0) for w G (—3, 1), Ri, is always a sink in the case A_|_, 
whereas we need to distinguish two subcases of the case A_, see Fig. 2 and the color-coding in 
Fig. 2. 

For (3 > f3\,(w) (i.e., in the + cases and the respective subcase of A_) there exists a fixed point in 
the interior of Su , which we call Cb : 
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C b The coordinates of C b are E + = 8+ l+^_ w) , M? = 36(i-^)[3g(i-gW(i + 3^)] 

This fixed point represents a solution that generalizes the Collins-Stewart solution (33); it reduces 
to (33) in the case /3 = 0. The point is a sink in S b ; this will be proved below by using global 
methods. When f3 converges to /3\,(w), the point C\, converges to IV Finally, the role of C|, in the 
type II state space X\\ is determined by 

[|:log S ] |Cb =-6E +|Cb <0, 

which implies that the fixed point C\, attracts orbits from the interior of X\\\ it is therefore repre- 
sented by a black dot. 

To complete the dynamical systems analysis of the set S b , we use the monotonicity principle with 
the function 

m? s n l ~ s , . 

z * - (T^F ' (58a) 



where 

C = 8+3/3(1— w) ' u ~~ (l+3ui)(l-c 2 ) 

This function is well-defined when j3 > — 1 (which is the difficult case); in fact, Z3 satisfies 



l+3w c _ 3(1— w)(c+f3)c \ 

- (i+3wxi-c^) ■ y^ ) 



Z' 3 = X Z 3 , where x-^3^y(S + - C ) 2 >0 (V/3 > -1) 

and x = if and only if S + = c. The implications of the monotonicity principle are analogous to 
those for the set <Sj. The results are summarized in the phase portraits depicted in Fig. 6. 

Global dynamics on the type II state space X\\. We are now in a position to complete the 
analysis of Bianchi type II models: By collecting the previous results concerning the four boundary 
components (52) of the type II state space X\\ and applying the monotonicity principle, we get a 
complete description of the dynamics in the interior of X\\ and thus of Bianchi type II models. 

Theorem 2. The qualitative behavior of typical (in fact, generic) orbits of the dynamical sys- 
tem (49) — which represents LRS Bianchi type II models — in the various anisotropy cases is the 
one sketched in Fig. 7. The past and the future attractors for the various cases are listed in Ta- 
ble 4- There exist non-generic solutions that are characterized by different asymptotic behavior; 
the complete list of possible a- and co-limit sets of orbits is given in Table 5. 

Proof. Consider the function Z4 : X\\ U Vn — > (0, 00) given by 



Z 4 = Mj 



-3/2 



1 - 2s 
which satisfies 

= -3 (£+ + |(1 + 3w)n) Z 4 and Z'l' = -\ZiM{ when (£+, 0) = (0, 0) . 

It follows that the function Z^ is strictly monotonically decreasing on AnUVn. By the monotonicity 
principle, the a-limit set 01(7) and the w-limit set ^(7) of every orbit 7 in the interior of X\\ must 
satisfy 

a( 7 ) cSjUTPi, w( 7 )c5 b . (59) 

It merely remains to study the flow induced on these boundary subsets to determine the possible 
a/w-limit sets. This is done by inspection of the phase portraits of Figs. 2, 5 and 6. 

The invariant structures on the boundary subsets (59) that are potential a/w-limit sets are fixed 
points, periodic orbits, and hcteroclinic cycles/networks. From Figs. 2, 5 and 6 we see that periodic 
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(a) B_, C_, D_ 



(b) A_ (-K/3</3 b ) 




Figure 6: Phase portraits of the flow on The comments of Fig. 5 apply analogously. 



orbits do not occur, so that the a/w-limit sets of orbits in X\\ must be fixed points or heteroclinic 
cycles/networks (if present). We begin by investigating the fixed points. 

Consider, for instance, the interior of X\, see Fig. 2. The only possible limit set in the interior of 
X\ is the fixed point F. In the — cases, F is a saddle in X\\ moreover, it is represented by a white 
circle in the subfigures of Fig. 2, hence it is a repellor in the orthogonal direction (i.e., into the 
interior of X\\). We infer that F is the a-limit set for a one-parameter set of orbits in X\\. which 
makes F a possible a-limit set, but not a generic repellor. In the + cases, F is a sink in X\\ since 
it is a repellor in the orthogonal direction, there exists exactly one orbit in X\\ that possesses F 
as an a-limit set. 

Inspection of the fixed points on dX\ yields only one (generic) source of orbits in X\\\ the Taub 
point TV This is in the — cases and in the case A_|_; see Figs. 2(a)-2(f), where T\, is a repellor 
within X\ and color-coded as a repellor in the orthogonal direction. The sinks we find are Qi, in 
the cases B_, C_, D_ and Ri, in the subcase — 1 < (3 < f3\,(w) of A_; in the remaining cases there 
does not exist any sink on X\. 

Turning to Sj we do not find any sources or sinks, see Fig. 5. On S\, there is the Taub point T[,, 
which is a source in the — cases and in the case A_|_, which is consistent with what we already 
know. The sinks are Qi, in the cases B_, C_, D_, R|, in the subcase — 1 < (3 < Pb(w) of A_, and 
C b in the subcase P^(w) < (3 < of A_ and in the + cases. 

While in the cases A-j-, B_, and C_, fixed points are the only possible a- and w-limit sets of 
interior orbits in Xjj, in the cases D_, B_|_, C_|_, and D_|_, there exist heteroclinic cycles on the 
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boundary which may belong to the past/future attractor. Let us study these cases in more detail. 

In the cases B_|_, C_|_, D_|_, the boundaries of the boundary subsets <Sj, X\, and <St, form a heteroclinic 
network, which is a set of entangled heteroclinic cycles, see Figs. 7(g) and 7(h). In particular, in 
the cases B_|_ and C+, the boundary of the set S$UXi, where the a-limit set of generic orbits must 
be located by (59), is a heteroclinic cycle, which is of the type (47); see the bold dashed lines in 
Fig. 7(g). This heteroclinic cycle remains as the only possible a-limit set of (generic) orbits in X\\. 
In the case D+, the set <Sj U dX\ consists of a heteroclinic network, cf. Fig. 7(h); the structure of 
this network is that of (47), where there exists a one-parameter family of possible paths from Tj 
to Qjj (since every orbit in 5j connects with Qjj, see Fig. 5(e)). It is not clear which part of the 
heteroclinic network acts as the actual a-limit set of orbits in Xu- In any case, it is immediate 
that the approach toward the singularity is oscillatory for generic solutions in Xu. 

To investigate the future asymptotics in the cases B_|_, C_|_, D_|_, we note that the boundary of 
S], forms a heteroclinic cycle connecting the fixed points T\, and Q|,. We claim that there are no 
interior orbits that have this heteroclinic cycle as the w-limit set, which leaves the fixed point C|, 
as the only possible aj-limit for all interior orbits. To prove this claim, consider again the function 
-Z3 given by (58); since c < 1 for (3 > 1, we may regard Z3 as a function that is defined on the 
entire state space Xu (except at £+ = c). Moreover, since S G (0, 1) for (3 > 1, Z3 vanishes on the 
boundary of S\>. Therefore, to prove our claim, it is enough to show that along any interior orbit, 
Z3 would be increasing if the orbit approached some point on the boundary of <S[,. In fact, by a 
direct computation one can verify that 



where 



rAlno-7 1 4(2+/3)(l+3u.)(£ + -c) 2 , fim 

Idr iu 6^3j | a=0 - 3c(l-w) (1-E+c) (fi+PiXP+Pa) ' ^ UU) 

Pi = W^T < 02 = • 
3(1 — w) " z 1— w 

Since /3i,/3 2 > 2, for w s (— 1,1), the r.h.s. of (60) is positive (unless E + = c), which establishes 
the claim. 

Finally, in the case D_ , there exists a family of heteroclinic cycles containing the fixed points 
Tjj and Qjj, see Fig. 7(a). The orbit that connects Tj with Qj is the semi-circle of dS$, whereas 
the orbits connecting Qj with Tj lie on the Bianchi type I boundary X\, cf. Fig. 2(a). To show 
that there do not exist orbits in Xu that have any of these heteroclinic cycles as the a-limit set 
a rather technical analysis is required; we merely outline the argument here. Assume that there 
exists an orbit 7 C Xu such that 0(7) = dS$ (or any other of the cycles). By assumption, there 
exists a sequence of intervals (t„ + i,t„), n G N, such that the solution 7(7") is in a neighborhood 
£+ G ( — 1, —1 + e) (and s > ^ — e') of Tj, i.e., —1 < S+(r) < — 1 + e for r n+ i < r < r n , n even, 
and in a neighborhood S + e (1 — e, 1) of Qj, i.e., 1 — e < < 1 for r„ + i < r < r„, n odd. In the 
limit r — > — 00, eq. (49b) decouples from (49a) and (49c), and we are able to compute the times 
At„ = T n — T n +i the solution spends in the neighborhood of Tj and Qj, respectively. We find that 
At„ is increasing in such a way that Ar n is strictly larger than the Ar„_i; more specifically, there 
exists c > 1 (which is independent of e) such that Ar„ > cAr n -i. Considering (49b) we see that 
\ — s increases (in the past direction) at least like e 6 *- 1-6 -* 1 " for r G (r ra +i,r„), n even, i.e., while 
the solution is in a neighborhood of Tjj. Conversely, i — s decreases (toward the past) at most 
like e~ 6T for r G {t„+i, t„), n odd, i.e., while the solution is in a neighborhood of Qj. Since, in 
the limit r — > — 00, the time the solution spends away from the fixed points becomes negligible as 
compared to the times Ar„, the behavior of s — ^ is represented by the alternating sequence 



^ — s ~ exp 



6((1 - e)Ar - An + (1 - e)Ar 2 - Ar 3 + (1 - e)Ar 4 



Since At„ > (1 — e) Ar„_i, n even, this behavior contradicts the assumption that 0(7) = 
since i — s -f* 0. This establishes the claims of the theorem. □ 
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Type 


Past Attractor 


Future Attractor 


B_,C_, D_ 


T b 




A_(/3< A) 


T b 


Rb 


A_(/3 b < 0), A , A + 


T b 


c b 


B+C+ 


Hetcroclinic cycle 


Cb 


D+ 


Heteroclinic cycle/network 


Cb 



Table 4: Past and future attractor of orbits in X\\, which represent LRS Bianchi type II cosmo- 
logical models. The exact solutions represented by the fixed points are given in Sec. 5 and in 
Appendix A. 



Interpretation of Theorem 2. In the + cases and in the subcase (3\, (w) < (3 < of A_ , the future 
asymptotic behavior of LRS Bianchi type II solutions is governed by the approach to the solution 
represented by the fixed point CV In Appendix A we compute this solution explicitly, see (103), 
and demonstrate that it constitutes a natural generalization of the Collins-Stewart perfect fluid 
solution (33). (Recall that for type II perfect fluid cosmologies, the Collins-Stewart solution is the 
future attractor.) In the subcase —1 < j3 < P\,(w) of A_, generic type II solutions are asymptotic 
to the solution represented by R[,, which is given as (101) in Appendix A; note that this is a 
non-vacuum solution because f2 7^ 0. Finally, in the cases B_, C_, D_, every generic solution is 
future asymptotic to the non-flat LRS Kasner solution (31), which is a vacuum solution. The past 
asymptotic behavior of anisotropic type II solutions resembles the behavior of perfect fluid solutions 
in the — cases and in the case A-|_: Towards the singularity, solutions behave asymptotically like 
the Taub solution. "Matter does not matter" in these cases, since the behavior coincides with that 
of vacuum solutions. However, in the cases B_|_, C_|_, D_|_, anisotropic matter matters. (Note that 
B_|_ is consistent with the energy conditions.) The approach to the singularity is oscillatory in 
these cases; the solution oscillates between the Taub family (30) and the non-flat LRS family (31). 
Since Q for these solutions, this is an example of asymptotic behavior that differs from the 
behavior of vacuum solutions. 

Remark. An example of a type B_|_ matter model is collisionless (Vlasov) matter for zero mass 
particles. In this special case, the existence of oscillations toward the past singularity for LRS 
Bianchi type II solutions was already observed in [28]. We refer to [8] for an extension of the 
analysis in [28]. 



9 Bianchi type IX: Setup 

9.1 The reduced dynamical system for Bianchi types IX (and VIII) 

For the Bianchi types VIII and IX, the structure constants h\ and 77,2 = 773 are non-zero; according 
to Table 1 we have h\ = 1 and 77,177,2 = =pl in type VIII and type IX, respectively. Setting h\ = 1 
and hifi2 = =Fl in the dynamical system (37) yields a system of five coupled equations; however, 
the constraints (40) imply that of the five variables in (37) only three are independent. 

Remark. We note that Hp 7^ 1 (i.e., H 7^ D), hence the normalization in (35a) is not the standard 
Hubble-normalization. For Bianchi type VIII (where 77477,2 = —1) the constraint (40b) implies 
that Hd > 1 (when we restrict our attention to expanding cosmological models). Bianchi type IX 
(where 77,1772 = +1), is characterized by Ho € (—1,1). 

It is not immediately evident, whether a three-dimensional dynamical system with suitable regu- 
larity properties can be constructed from the constrained system (37). To achieve this aim, in a 
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(a) D_ 



Qb 

(b) B_, C_ 




Figure 7: Phase portraits of Bianchi type II orbits. Bold lines depict the behavior of typical 
orbits. Dashed lines are orbits which lie on the boundary of the state space and form heteroclinic 
cycles/networks. In the cases B_|_, C+, D+, the heteroclinic cycle/network depicted with bold 
dashed lines is the past attractor of interior orbits. In all other cases, the past attractor is the 
fixed point TV The future attractor is a fixed point in all cases. There exist non-generic orbits 
whose behavior toward the past or toward the future is different; see Table 5. 
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Fixed point 




Dimension of unstable [stable] manifold intersected with X\\ 


D_, C_ 


B_ 


A_ (/3 < A) 


A_ (/3 b < p) 


A+(/3</3») 


A+(/3, < /3) 


B+, C+, D+ 




3 




3 


3 


3 


3 




Qb 


[3] 














R-b 






[3] 




2 


2 




Rj 










2 






F 


2 




2 


2 


1 


1 


1 


c b 








[3] 


[3] 


[3] 


[3] 




1 




1 


1 


1 






het. cycle 














3 



Table 5: Possible a- and w-limit sets of LRS Bianchi type II cosmological models. Each Bianchi 
type II solution corresponds to an orbit in X\\\ hence the intersection of the unstable [stable] 
manifold of a fixed point with X\\ yields the set of type II solutions converging to that fixed point 
as t —t [t — > oo]. If this set is three-dimensional, the fixed point is a source [sink]; if the dimension 
is two, there exists a one-parameter set of orbits converging to the fixed point; if the dimension 
is one, there is only one orbit with that property Center manifold analysis shows that the cases 
A -|- and A _ behave like A_|_ (0 < j3 < /Jj) and A_ (/% < j3 < 0), respectively. 



first step, we closely follow [16]: We use (40b) to express M 2 in terms of Hp and Mi, so that (40a) 
becomes 

3fMi a (l-.H£) r i_ fry , 3|l-ffl,| r i_ 



Ml J V M\ 

which expresses s in terms of Hjj and M\ . The independent variables we choose are thus Hp , E + , 
and Mi , and the dynamical system (37) reduces to 

H' D = -(I - H 2 D )(q - H D Z + ) , (62a) 
£'+ = -(2 - q)H D ^+ - (1 - H 2 D ){1 - E 2 + ) + | h\M\ + m (u(s) - w) , (62b) 
M[ =Mi(qH D -AV + + {l-Hl)H + ) , (62c) 

where s is regarded as a function of Ho and Mi, see (61). The Hamiltonian constraint (39'), 

n = 1 - Y? + - ± nfM'l , (63) 

is used to solve for Q, and q is given by q = 2Y? + + |(1 + 3w)Jl. Note that the system (62) is 
identical for Bianchi type VIII and type IX; however, the state spaces are different: Hq > 1 for 
type VIII and H D £ (-1, 1) for type IX. 

Anticipating the analysis to follow, let us note that the system (62) does not have the regularity 
properties that are necessary to perform a satisfactory dynamical systems analysis. This is due to 
the fact that the variable s and thus the function u(s) in (62b) does not have a limit when H D — > 1 
and Mi — > 0. To remedy this problem, it will become inevitable to introduce yet another system 
of coordinates in a neighborhood of the 'bad' part of the boundary. 

Remark. Let us comment on the role of the boundaries. In the case of Bianchi spacetimes where the 
matter is assumed to be vacuum or a perfect fluid, the Bianchi types VIII and IX are characterized 
by the maximal number of 'true degrees of freedom': The dynamical system describing the dy- 
namics of (diagonal) vacuum [perfect fluid] models of Bianchi type VIII or IX is of dimension four 
[five]; in the LRS case the dimension is two [three]. In contrast, for the lower Bianchi types, one 
or several degrees of freedom are 'gauge', and the dynamical systems describing the lower Bianchi 
types are of lower dimension. In particular, one finds that the state spaces and the dynamics of the 
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lower Bianchi types appear on the boundary of the state spaces of the higher Bianchi types. This 
hierarchy of invariant subsets (which is called the Lie contraction hierarchy, see Appendix B.3) 
plays a fundamental role in the dynamical systems analysis of cosmological models, see [14]. In the 
case of Bianchi spacetimes where the matter is assumed to be anisotropic, such a hierarchy exists 
only partially. While the Bianchi type I state space appears as the boundary subset of the type II 
state space, see Sec. 8, the type II state space is not directly related to a boundary subset of the 
type VIII/IX state space. This is a simple consequence of the fact that the two state spaces have 
the same dimension. 

9.2 Kantowski-Sachs and LRS Bianchi type III models 

In this paper we refrain from giving an analysis of the dynamics Kantowski-Sachs and LRS Bianchi 
type III models. However, Kantowski-Sachs and LRS Bianchi type III play a certain role in the 
analysis of the flow on the boundary of the state spaces of Bianchi type VIII and IX. 

As observed in Appendix B, the spatial Ricci curvature of Kantowski-Sachs and LRS Bianchi 
type III models are obtained from the Ricci tensor of Bianchi class A models by formally setting 
n\ to zero and n\fii to +1 or — 1, respectively. As a consequence, the same is true for every system 
of equations describing these models. In particular, the (reduced) dynamical system representing 
Kantowski-Sachs models and LRS Bianchi type III models is obtained from (37) by formally setting 
n\ to zero and n\fi2 to +1 or —1, respectively. In complete analogy to the considerations in the 
type VIII and the type IX case, we are led to the dynamical system (62) and the constraint (39'), 
where 

nj = 0. 

Again, the variable s is given in terms of Mi and Hp by (61). Note that the system (62) is identical 
for Kantowski-Sachs and LRS type III; however, the state spaces are different: 

*ks = {(H D , S +> Aft) | H D G (-1, 1) , Aft > , Y? + < l} , (64a) 
X m = {{H D , £+, Mi) \H D >1, Ah > , < l} . (64b) 

Remark. In the case where the matter is a perfect fluid, i.e., u(s) = w, the equation for Mi 
decouples from (62) and the reduced dynamical system for Kantowski-Sachs models and LRS 
type III models is merely two-dimensional, the equations being (62a) and (62b). However, in the 
case of an anisotropic matter source, the system (62) is a coupled system, since the variable s, 
which appears in (62b), depends on Mi by (61). 

10 Bianchi type IX: Analysis 

The dynamics of LRS models of Bianchi type IX are represented by the reduced dynamical sys- 
tem (62), where Hrj 6 (—1,1) (and n\ = 1). Accordingly, the state space of LRS Bianchi type IX 
models is 

Xix = {(H D , £+, Mi) | H D G (-1, 1) , Mi > , Y? + + ^M 2 < l} , 

which is relatively compact. The state space Xjx is depicted in Fig. 8; like Xu it has the form of 
a tent; note, however, the difference in the set of variables. 

The state space Aix admits a natural decomposition into two 'halves' that are separated by a 
plane, i.e., 

Xix = XX U X® x U X lx , 
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Figure 8: The state space X\ X . 



where 

X+ = {H D > 0} c X 1X , = {H D = 0} c X IX , X 1X = {H D < 0} c X IX . 

The plane X® x acts as a 'semipermeable membrane' for the flow of (62), since 

H' D | ffD =o = -(2S 2 + + |(1 + 3w)n) < . 

This a subset that is past invariant under the flow of (62), while X^ is future-invariant. 
The physical interpretation is straightforward: A type IX model that is expanding at t = to, i.e., 
H(to) > (<=> Ho(to) > 0), must have been expanding up to that time, i.e., for < t < t ; 
conversely, a model that is contracting at t = ti(Le., H(t\) < 0) must continue to contract Vi > t\ 
(which eventually leads to a big crunch); finally, a model that satisfies H = at some time, is a 
model that starts from an initial singularity, expands to a state of maximum expansion (H = 0) , 
from which it then recontracts to a final singularity. (This "closed universe recollapse" behavior is 
the well-known behavior of vacuum and perfect fluid solutions of Bianchi type IX, sec [34] and [16]. 
However, there exist anisotropic matter models that deviate from "closed universe recollapse" , 
see [7] and the results presented in Sec. 11.) 

Remark. The dynamical system (62) is invariant under the discrete symmetry 

H D — » —Hd , S + -> -S + , t — > -t ; (65) 

the flow in X^ is the image of the flow in X^ under this discrete symmetry. In particular, it suffices 
to analyze the asymptotic behavior of solutions in X^; the asymptotic behavior of solutions in X^, 
follows by applying the discrete symmetry (where the roles of a- and w-limit sets are reversed). 

10.1 The boundaries of A'ix 

The boundary of the type IX state space X\ X consists of the invariant sets 

<S B = {H D = 1 ; -1< S + < 1, < Mi < ^12(1 - Y? + ) } , (66a) 
Six = {Mi = ; -1< £+ < 1, -1< H D < 1 } , (66b) 
Vix = {^ = ; -1< £+ < 1, -1< H D < 1, Mi = ^12(1 - } . (66c) 
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The symbol Sj is used with slight abuse of notation; we refer to the discussion of (69) below 
(and to the analogous comment following (52)). There exists a second copy of the set (66a) with 
Hd = — 1; since this set is merely the image of S% under the discrete symmetry (65), we may 
restrict our attention to S$. Analogously, concerning the 'base' 23ix, it is enough to restrict our 
attention to the past invariant half = Bix H {Hd > 0}, see Fig. 8. 

The dynamical system (62) admits a regular extension to each of the invariant boundary sub- 
sets (66). However, the system can not be extended to the entire boundary dXjx] this is because 
the variable s and thus the r.h.s. of (62b) does not possess a well-defined limit when (Hd, E+, Mi) 
converges to a point on the closure of the line 

Cj = dS t n 8B 1X ; (67) 

the use of the subscript I will become clear in Subscc. 10.2, when we remedy this defect of the 
system (62). 

The first step of our analysis is to study the flow that is induced on the boundaries (66) of X\x, 
where the dynamical system (62) admits a regular extension. 

The vacuum boundary Vix- The dynamical system (62) admits a smooth extension onto Vrx 
which is given by setting = in (62). Since — and thus p = on Vrx, the orbits on this set 
represent type IX vacuum models. The induced system is 

H' D = (l-H%)(H D -2S+)Ji + , (68a) 

E+ = 0- ~ S +) t 2 + (! " S +) + ( H d - £+) 2 ] 5 (68b) 

the variable Mi is given by Mi = 12(1 — £+). The closure of the vacuum boundary is a square, 

Vix = [— 1, 1] x [—1,1]. Eq. (68b) entails that £+ is strictly monotonically increasing whenever 
E + 7^ ±1. It follows straightforwardly that the vertices of the square are the only fixed points; two 
of these fixed point satisfy Hd = 1: 

T The Taub point T is given by (H D , £+) = (1, -1). 

Q The non-fiat LRS point Q is given by (H D , E + ) = (1, 1). 

The Taub point T corresponds to the Taub solution (30) and Q to the non-flat LRS Kasner 
solution (31). The reason for the lack of a subscript j or ^ will become clear in Subsec. 10.2. (The 
fixed points that correspond to T and Q on the subset Hd = — 1 are (i?_o,E+) = (—1,1) and 
(Hd, E+) = (—1, —1); the former is a Taub point, the latter a non-flat LRS point. The properties 
of these fixed points follows from the discrete symmetry (65).) 

Since £+ is a strictly monotonically increasing function on [—1,1] x (—1,1), the monotonicity 
principle applies. In combination with the fact that H' D < on the boundaries Hd = ±1, we 
arrive at the conclusion that the Taub point T is the a-limit set of every orbit in Vix (while the 
Taub point (—1, 1) is the w-limit set); see Fig. 9. 

The side S$. The dynamical system that is induced by (62) on the side <S(j of the type IX state 
space is identical to the system (53): 

EV = \Ml{2 - E + ) - §0(1 - «,)(£+ - |) , (69a) 
M[ = Mi [|(1 - w)Y? + - 4E+ + |(1 + 3w)(l - j^M?)] , (69b) 

where = 1 — E^_ — j^Mf . The system (69) on iSj possesses a smooth extension to <Sj (which does 
not contradict the fact that (62) cannot be extended to <Sj). 
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Figure 9: The flow induced on Vrx- 



The phase portraits for the various cases are depicted in Fig. 5. Note that the fixed point Cjj is 
represented by a white dot in the figures, which means that it acts as a repellor in a direction 
transversal to Sj, i.e., into the interior of X\x- To see this we use (62a) and find 

[£ log |f - H D \] lCf = 2(2E^ - E+ + \ (1 + 3w)Si) | c „ = nff^gg^ > ° ■ (™) 

The vertices of the boundary dS$ coincide with the fixed points Tj and Qj, see Fig. 5; in addition, 
in the cases A, B, C, the fixed point Rtj is present. The significance of these fixed points in the 
type IX context will become apparent in Subsec. 10.2. 

Remark. The solutions of (69) on the side S$ are type II solutions associated with an anisotropic 
matter source with 'rigid' rescalcd principal pressures, i.e., u(s) = const. For such an anisotropic 
matter source, cq. (49b) decouples from the system (49); the remaining equations (49a) and (49c) 
coincide with (69). Analogously, the orbits of (71) on the base Byx, to be discussed next, represent 
solution of Kantowski-Sachs type, see (116), where again the rescaled principal pressures are 'rigid'. 
For such an anisotropic matter source, eq. (62c) decouples from the system (62). (Recall that h\ is 
set to zero for the Kantowski-Sachs system.) Therefore, the associated Kantowski-Sachs models are 
described by the reduced system consisting of (62a) and (62b), which is identical to the system (71) 
induced on B\x- 

The base Bix- The dynamical system (62) extends smoothly to the base Bjx of the state space 
(but not to £>ix)- The induced system is obtained by setting Mi = in (62): 

H' D = -(1 - H%) [2 - |(1 - w)(l - E 2 + ) - H D E + ] , (71a) 

£'+ = -(1 - £+) [(1 - Hi) + |(1 - w){H D H + + /?)] . (71b) 

The induced system (71) on B\x possesses a smooth extension to Byx = [—1, 1] X [—1, 1] (which 
does not contradict the fact that (62) cannot be extended to B\x)- 

The boundary of B\x contains several fixed points. 

T b The Taub point T b is given by (Hd, £+) = (1, -1). 

The non-flat LRS point Qt is given by (Hd, £+) = (1,1)- 
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Remark. The subscripts (i.e., \, or j) are chosen with care; their significance in the present context 
will become clear in Subsec. 10.2. Although the points T on dVix, Tj on dS$, and T b on dB\x 
have formally the same coordinates (i.e., Ho = 1, E + = —1, Mi = 0), these points cannot be 
identified. We will see in Subsec. 10.2 that the Taub points T b , Tjj, T provide different (but 
equivalent) representations of the Taub solution (30). The reason for this is that the system (62) 
cannot be extended to the line C\. While the extension of the system (71) on B\x to the boundary 
subset Hd = 1 yields eq. (46^), the extension of the system (69) on S§ to (the base of) <9Sj leads 
to eq. (46j). Unless (3 = 0, these two equations are different. 

If and only if \(3\ < 1 (i.e., in the cases A_|_ and A_), the boundary of B\x contains the point R^. 



Rb This fixed point has the coordinates (Ho, £+) = (1, — /3). 



Applying the discrete symmetry (65) we obtain an analogous fixed point with Ho = — 1. In 
Appendix A we calculate the explicit solution associated with Rj, and give its interpretation. 

Searching for interior fixed points in 2?ix we find a fixed point P in B^ (and an analogous point 
in generated by the discrete symmetry). 

2+3,3(1-™) v _ l+3w 



P The coordinates of P are H D = ^w-^; E 



^/(l-3ui) 2 +6/3(l-to) ' + ^/(1-3™) 2 +6/3(1-kj) ' 

There are, however, two requirements for the fixed point P to exist: First, w must satisfy 

-§ < w < « -0.244 (72a) 

and, second, (3 must be in the range 



/3_</3</3 + with P ± = 1 % 3 1 +l) J^l. (72b) 



3(1— w) 

The conditions (72) ensure that Hp | P and S + | P are real (and positive) numbers strictly less then 
1. The range of admissible values of w and (3 is depicted in Fig. 10; since — | < /3_ < — 2 + \/3 
and —2 + < (3+ < 0, cf. Fig. 10, we find that the conditions (72) distinguish a subcase of A_, 
which we call the Bianchi type IX special case and denote by A : i^. 

An alternative representation of (72b) in terms of w_, see (19), is 

i (l + 6w - v/-3 + (1 - 3w) 2 ) < v- < \ (l + 6w + v / -3+(l-3u>) 2 ) , (72b') 

cf. Fig. 11. Finally, note that P coincides with R|, in the limiting cases (3 = (3±. The exact solution 
corresponding to the fixed point P is of Kantowski-Sachs type and is derived in Appendix A. 

The local dynamical systems analysis of the fixed points is straightforward. The role of the fixed 
points Tt, and is determined by the flow on the boundaries of Six, which is obtained in a 
straightforward way; we simply refer to Fig. 13. The fixed point Rt, is an attractor on the line 
Hd = 1; in the transversal direction we have 

[§p log (1 - H D )] |Rb = (1 + 3w) + 2/3 + 3(1 - w)f3 2 , (73) 

which is negative if and only if w satisfies (72a) and (3 satisfies (72b). Hence, in the Bianchi 
type IX special case A 1 ^, the fixed point Rj, is a sink in Bix', in all other cases, R[, is a saddle. (In 
the marginal case (3 = (3±, this does not follow from the local analysis since R[, is not hyperbolic 
in these cases. However, the statement is still correct, as will be proved by using global methods.) 
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Figure 10: The Bianchi type IX special case A 1 * is defined by conditions on w and the anisotropy 
parameter (3. The admissible values are depicted as the gray region. The axes of the diagram are 
(3 = and w = -0.23. 



The fixed point P which exists in the Bianchi type IX special case A 1 ^ is a saddle in Bix- This is a 
direct consequence of the linearization of the dynamical system (71) at P; in fact the determinant 
of the linearized system is 



27(l-w) 3 



(1-3w) 2 +6I3(1-w) 



which is negative because the three terms in brackets are positive under the conditions (72a) 
and (72b) of case A™. Accordingly, the linearization is indefinite and the point P is a saddle. 

Since the fixed point P lies within Bix (and not on the boundary Ho = 1 like T^, Qt,, Rt,), the 
dynamical system (62) on X\x extends smoothly to a neighborhood of P in Xyx- Using (62c) we 
find 



[£l0gM!] p = -- 



3(l+3tu) 



< 0. 



(74) 



y/(l-3w) 2 +eP(l-w) 

Therefore. P is a saddle in X\x with a two-dimensional stable manifold and a one-dimensional 
unstable manifold. We reemphasize that analogous statements for the fixed points T[,, Q^, R], are 
unavailable in the present dynamical systems formulation, since the system (62) does not extend 
to the line C\, cf. (67). In Subsec. 10.2 we will remedy this defect. 

Let us turn to the global dynamics on 2?ix- 

Lemma 1. The phase portraits of the dynamical system (71) on the base Bix for the various 
anisotropy cases are as depicted in Fig. 13. 



Proof. The eq. (71a) for Ho can be written as 



H' D = -(1 - Hi) |(1 + 3w) - H D ^+ + |(1 - w)T^ 



(75) 



For Ho > and E + < 0, the expression in brackets is manifestly positive. For positive S + the 
expression attains its (unique) minimum for the value Ho = 1; in that case we obtain a second 
order polynomial in (— £+), which coincides with (73) (modulo a factor of 2) and is therefore non- 
negative if w > 1 ~ 3 V ^ . We conclude that the expression in brackets is positive for all Ho G [0, 1), 
£+ G (-1,1), if w > 1=&. 
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Figure 11: The admissible values of w and V-/v+ for the Bianchi type IX special case A 1 *. The 
axes are v± = and w = —0.23. From this figure, in combination with the considerations of Sec. 4, 
it is evident that the dominant energy condition and the strong energy condition are satisfied, since 
\v±\ < 1 and w > -1/3. 



This implies that, in the case w > — , Ho is strictly monotonically decreasing on the sub- 
set Ho > of Bix, i.e., on B^ x . Applying the symmetry (65) we see that Ho is then strictly 
monotonically decreasing on the entire state space Byx- Taking into account the simple structure 
of the flow on dBix, we conclude that the a-limit set of an orbit in Six must be one of the fixed 
points with Ho = 1: T|>, Qi,, or R|, (which requires |/?| < 1, i.e., one of the A cases); the former 
two possibilities represent the past asymptotic behavior of typical orbits, which is because Rt, is 
a saddle when w > 1- g (except in the case w = x ~g , (3 = (3±, where it is a center saddle, as 
follows from a center manifold analysis). By applying the symmetry (65) we obtain an analogous 
statement for the possible w-limit sets. In the A cases, when the fixed point Rt, exists, there is 
one special orbit emanating from Rt, that acts as a separatrix in the state space B\x- To obtain 
the w-limit point of this special orbit we consider the function r\, = (3 Ho + S + on Byx] the set 
r\, = represents a straight line in Byx that connects the fixed point Rt, with its counterpart on 
Ho = — 1. A straightforward computation shows that 

W)|, b= o = -(1-^)(1 + 2/3), (76) 

which is negative for (3 > —1/2, zero for j3 = —1/2, and positive for j3 < —1/2. Therefore, the 
subset r\, < of B\x is future- invariant if (3 > —1/2, while r\, > is future-invariant if j3 < —1/2. 
(The set r\, = is an actual orbit in the case (3 = —1/2). The local analysis of the system (71) at 
the fixed point Rt shows that the unstable manifold of Rt, and thus the separatrix orbit emanating 
from R[, is contained in the future- invariant subset. From this observation it is straightforward to 
determine the w-limit point of this orbit. We conclude that the structure of the flow is indeed as 
depicted in Fig. 13. 

In the case — ^ < w < 1 ~ 3 V ^* . the monotonicity of Ho on B\x is violated; however, the property 
H' D > is restricted to a rather small region of B\x, sec Fig. 12. We distinguish two subcases: 
(3 £ (/3_, (3 + ) and (3 £ (/3_, j3+), where the latter coincides with the special case A™. First, consider 
j3 £ (/3_,/3 + ). In this case, there does not exist a fixed point in the interior of Byx- This fact 
excludes the presence of periodic orbits and heteroclinic cycles in B\x- (For a two-dimensional 
system, a periodic orbit necessarily encircles a fixed point which is a sink or a source.) The 
absence of these structures in Byx and the Poincare-Bendixson theorem (see Appendix C) imply 
that an orbit in B\x must leave every compact subset of Six; in other words, the a- and w-limit 
set must intersect the boundary. Taking into account the flow on the boundary it follows that the 
possible a-limit sets are the fixed points on Ho = 1; the possible w-limits sets are the fixed points 
on Ho = —1. Using again (76) we conclude that Fig. 13 depicts the correct qualitative behavior 
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Figure 12: These figures depict a neighborhood of the boundary Hp = 1 on Six- On Hr> = 1 
(which is the vertical line) there are the fixed points T[,, Qt,, and, in the A cases, the fixed point 



R b . If 



< w < 



l-\/3 
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there exists a small region in Six, where H' D > 0. If (3 $ (/3_,/3+), which 



is the case for the figure on the l.h. side, this does not have any influence on the qualitative global 
dynamics on Six- However, if /3 € (/?_, f3+), which corresponds to the special case A 1 *, then there 
exists a fixed point, P, in Six, and R[, is a sink. It is not difficult to construct a future-invariant 
set that is a neighborhood of R[, and contains an orbit P — > R b ; see the figure on the r.h. side. 



of the flow on Six (where, however, Hd is not necessarily monotonically decreasing). 

Second, consider (3 £ (/3_, /3+), i.e., the Bianchi type IX special case A 1 —, in which the presence of 
the fixed point P slightly complicates matters. However, the fixed point P (as well as its counterpart 
in the subset Hd < 0) is a saddle; therefore, there do not exist any periodic orbits in Six and 
the following argument excludes the existence of hcteroclinic cycles as well. Since P lies on the 
boundary of the region H' D > 0, it is not difficult to prove that one of the two orbits emanating from 
P must converge to R b . (The first step is to construct a future-invariant subset of the set H' D > 0, 
whose closure contains the sink R b . Taking into account the flow on the boundary Hd = 1, we 
find that R[, is the attractor for every orbit in this future-invariant set; see Fig. 12. The second 
step is to show that one half of the unstable subspace of P is contained in this future- invariant set. 
We conclude that one orbit emanating from P is contained in the future-invariant set and thus 
converges to Rt, which proves the claims.) The w-limit set of the second orbit emanating from P 
must be contained in the future-invariant half Ho < of the state space Six- More specifically, this 
orbit converges to the fixed point (—1, —1) (which represents a non-flat LRS point) on Hu = — 1. 
(To prove this claim we note that r^p < 0, i.e., P is contained in the subset r b < of the state 
space Six; this set is future-invariant since (3 > (3- > —1/2. The only possible attractor in the 
future- invariant set {r\, < 0} set is the fixed point (—1,-1).) Collecting the results we have proved 
that Fig. 13(e) depicts the correct qualitative behavior of the flow on Six in the special case A 1 —. 
This complete the proof of the lemma. □ 
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(a)B_,C_,D_ (b) A_ (p < -1/2) 




(g) B + , C+, D + (h) Coordinates 



Figure 13: Phase portraits of the dynamical system (71) on Z?rx- Typical orbits are represented 
by continuous lines, non-generic orbits by dashed lines. The color-coding of the fixed points T\,, 
Qb, Rb makes sense only in the context of the dynamical systems formulation of Sec. 10.2; the 
coordinate r of Subfig. (h) is used in that section. 
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(a) The space 3^ (b) A schematic depiction of f 3^ix 

Figure 14: The space 3^ T ~t and its boundary. 



10.2 The Bianchi type I boundary 

In this subsection we introduce the tools that are needed to study the behavior of solutions of the 
dynamical system (62) in a neighborhood of the line C\ = dS§ n dBix, cf. (67). A careful analysis 
is necessary since the system (62) does not admit a smooth extension to that part of the boundary 
of the state space Abe- 

To remedy this defect of the system (62) wc introduce a set of 'polar coordinates' centered on the 
line L\. Let 

Af? = 3rsintf, (77a) 
l- J fff, = 2rcostf , (77b) 
£+ = unchanged . (77c) 

The transformation of variables from (i?D,Mi,£ + ) to (r, where r > and < •& < ^, is 

a diffcomorphism on the domain (Hd,Mi, £+) £ XX. We define the domain J^x °f the variables 
(r, £ + ) to be the preimage of the state space XX^ under the transformation (77). We obtain 
rcostf < 5 from (77b); the constraint 1 - £+ - i Mf = £1 > implies rsini? < 4(1 - £^_). 
Therefore, yX^ can be written as 

y& = { r > 0,1? G(0, |), £ + G(-l ) l)|r<min[ f ^ ) ^^] }, (78) 

see Fig. 14(a). By construction, the state space J^x is P as t invariant under the flow of (79) (because 
XX^ is past invariant). 

In the new variables, the dynamical system (62) takes the form 

r' = 2r (H D (q - H D T,+) - 3£+ sin 2 §) , (79a) 
= -3£+ sin(2tf) , (79b) 
£' + =rsini?- 1 + (H D — £+) 2 + Hr>Y, + {q — Hd^+) + 3£1 (i*(s) — w) , (79c) 

where q = 2£ 2 f _ + i(l + 3u>)f2 and fi = 1 — £l — \ r sini?. In addition, ifo and s are regarded as 
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functions of r and d in (79), 



1 tan t9 

H D = VT=*FZZ*, s =2lTw- (80) 

Note that these are smooth functions of r and "& on the state space 3^y- In particular, 

(9s 1 11 1 (9s 1 

= 2s(l-2s) = ■ — 777777 ; hence - < ^ < - Vi? G [0, f ] (81) 

d?9 sin(2#) v 7 2 1 + sin(2i?) 4 ~ <9tf ~ 2 1 ' 2 J y ' 

and <9s/<9i? = | at t? = and i9 = |. 

In contrast to the system (62), the dynamical system (79) on the state space y^ admits a regular 
extension to the boundaries (where we assume Ho > const > 0). The reason for this is that s and 
thus the function u(s) in (79c) is well-defined for each point on the boundary of the state space. 

The side S$. In terms of the variables (77), the set S$, cf. (66a), corresponds to 

Sf = {tf = | ; -K £+ < 1, < r < 4(1 - S 2 ,) } . 

(We use the symbol 5h with slight abuse of notation.) The induced dynamical system reduces 
to the system (69), where Mi is replaced by y/3r. The phase portraits for the various cases are 
depicted in Fig. 5. 

The base Bix- In terms of the variables (77), the (past invariant half of the) set Z?rx, cf. (66b), 
corresponds to 

#ix = = o ; -i < s + < 1. < r < \ } • 

It is an invariant subset of the system (79) on J^ix as well. The induced system on Bf x is equivalent 
to the system (71) on the past invariant half of Six, when we set Hd = y/l — 2r. The phase 
portraits for the various cases are depicted in Fig. 13, where we restrict our attention to the past 
invariant (i.e., the right) half. 

The vacuum boundary. Setting SI = in (79) yields the (past invariant half of the) vacuum 

subset. The vacuum subset of y^ x is different from Vj'x, which is depicted in Fig. 9. The main 
difference concerns the fixed points: While there are two fixed points, T and Q, on the past 
invariant half of the vacuum subset V IX of Af IX , there are four fixed points on the vacuum subset 

T b The Taub point T b is given by (r, •&, £+) = (0, 0, -1). 

Tjj The Taub point T, is given by (r, S+) = (0, f , -1). 

Q b The non-flat LRS point Q b is given by (r, tf, E+) = (0, 0, 1). 

Qt) The non-flat LRS point Q B is given by (r, ??, S+) = (0, f , 1). 

The transformation (77) duplicates the fixed points T and Q; more specifically, there exist two 
orbits on the vacuum subset of J 7 do an orbit T b — > Tj and an orbit Qt, <— Qj, that are collapsed to 
the fixed points T and Q on V^ x . This becomes apparent when we consider the fourth boundary 
subset of y^x, which does not have a direct correspondence on X IX . 

The Bianchi type I boundary X \. There exists an additional two-dimensional boundary subset 
of y^: The set r = 0. 



Xi = {r = ; -1< £+ < 1, < ■& < f } 
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This boundary subset can be regarded as the preimage of the line C\ under the map (77). Since 
the map (77) is not injective for r = (because (r,??) = (0, i9) i— » (Mi, Hp) = (0, 1), W), the set 
X\ is two-dimensional; the introduction of the coordinates (r, $) generates a blow-up of the line 
Cj. The dynamical system (79) extends regularly to the subset X\; the extension reads 

•& = -3S+ sin(2tf) (82a) 
EV = -3(1 - ) [|(1 - - («(«) - to)] . (82b) 

The variable s is given in terms of i? by (80); from (81) it is immediate that the system (82) 
coincides with the Bianchi type I system (42); therefore, we call X\ the Bianchi type I subset. 
The phase portraits for the various cases are depicted in Fig. 15; the difference between Fig. 2 and 
Fig. 15 is in the color-coding of the fixed points G> and R|,. Note in particular that attracts 
orbits from the orthogonal direction in the Bianchi type IX special case A™, see Fig. 15(c); this 
follows by computing 

f T [logr] |Rb = 3(1 - w)p 2 + 2/3 + (1 + 3w) , 
which is negative in the case A 1 * and thus consistent with Fig. 13(e). 



11 Bianchi type IX: Results 

In this section we use the building blocks of Sec. 10 to derive the main results on the global dynamics 
of LRS Bianchi type IX models. The physical interpretation of these results is stated explicitly. 

Let 7 denote an orbit of the dynamical system (62) on Xix, which represents an LRS Bianchi 
type IX model. There exist two fundamentally different cases in the context of Bianchi type IX 
dynamics. 

• The recollapsing case: 7 n X^ x ^ (=> 7 (jt X^). 

• The (forever) expanding case: 7 C X^ (7 n X® x = 0). 

In the recollapsing case there exists To > such that 7(7-) £ X^ (i.e., Hd > 0) Vr e (— 00, to), 
7(t ) G (i.e., Hd{tq) = 0), and 7(7-) G X^ (i.e., Hd < 0) Vr < r < 00. The cosmological 
model expands from an initial singularity, reaches a point of maximum size, and then recollapscs 
to finish in a big crunch. 

In the expanding case (non-recollapsing case) we have 7 C X^ i.e., 7(7") G X^ (and thus 
Hd > 0) Vr G (—00,00). In this case the cosmological model expands forever. 

A third case, the collapsing case, arises from the expanding case by the discrete symmetry (65) 
and corresponds to forever contracting models. 

Lemma 2. Let 7 C X\x be a recollapsing or expanding orbit of the dynamical system (62). Then 

a(i) C So • 

In the recollapsing case, ^(7) is contained on the image of Sj under the discrete symmetry (65), 
in particular, Hd — ■* — 1 as r — > 00. In the expanding case, to(j) C Bix- 

Remark. Note that both Sj and Bix include the type I set C\. 
Proof. Consider the C 1 function Z5 : X\x U Vix - ¥ K given by 

A/ l/3 

= H D p/3 • (83) 
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Figure 15: Phase portraits of orbits on the Bianchi type I subset X\ of D^rx- The fixed points are 
color-coded: A black fixed point is an attractor, a white fixed point is a repellor in the orthogonal 
direction. This figure differs from Fig. 2 in the color-coding of the fixed points Qt> an d Rb- 



Since Z5 satisfies 

z 5 = (2^ + |(1 + 3w)n) 

and 

^5" = ~ ^_%v* V + ) 2 if fi) = (0, 0) , 

this function is strictly monotonically decreasing on AixU Vrx- The monotonicity principle implies 
that the a- and the w-limit set of every orbit in Xjx must be contained on 9A4x\Vtx (which is the 
union of *8ix, <5fl, and the counterpart of 5j with Hd = — 1). 

For a recollapsing or an expanding orbit we have Z5 > for some time t_ . Because Z§ is zero on 
Six, we infer that 0(7) n £?ix = 0, The first claim of the lemma ensues. 

For a recollapsing orbit we have Z5 < for some time t+. Therefore, 07(7) n Six = 0, because Z5 
is zero on Six- For an expanding orbit, Z§ is positive for all times; this leaves 0^(7) C Bix as the 
only possible asymptotic behavior. This concludes the proof of the lemma. □ 
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To proceed we make use of the dynamical systems formulation of Subscc. 10.2. An orbit 7 of the 
dynamical system (79) on y^ represents the expanding phase of an LRS Bianchi type IX model; 
this is because the past invariant subset X^ of X\x (which is defined by Hp > 0) and 3^rx are 
diffcomorphic. (In the expanding case, there is but the expanding phase; hence, in that case, 7 
represents the entire type IX solution.) Since the dynamical system (79) on y^ possesses a regular 
extension to the boundaries of the state space, we are able to apply standard methods from the 
theory of dynamical systems. 

Lemma 2'. In the context of the formulation of Subsec. 10.2, the results of Lemma 2 imply that 
C «S» U Afj. Furthermore, a; (7) is either contained in Bix U X\ (expanding models) or there 
is a big crunch (recollapsing models). 

Theorem 3. In the Bianchi type IX special case A 1 * defined by (72), recollapsing solutions and 
forever expanding solutions occur on an equal footing. The qualitative behavior of typical solutions 
is the one sketched in Fig. 17. In all other anisotropy cases, every model is a recollapsing model. 

Proof. Assume an anisotropy case different from A 1 *. If there were an expanding model 7, then 
^(7) C Bix U X\ by Lemma 2'. Inspection of the flow on Bix reveals that there do not exist any 
potential w-limit sets on Bix', in particular, the stable manifold of every fixed point has an empty 
intersection with the state space J^ixi see Fig- 13. The same is true for X \, see Fig. 15. The only 
slightly non-trivial case is the case D+, where dX\ is a heteroclinc cycle. However, we may use the 
function Z\, see (48): If a; (7) = dX\, then Z\ — > 00 along 7 as r — * 00. However, a straightforward 
computation shows that Z\ is decreasing in the neighborhood of the fixed points of dX\ which 
contradicts the divergence of Z\. 

Second, assume the Bianchi type IX special case A 1 *. Figs. 13(e) and 15(c) imply that there are two 
potential w-limit sets for orbits 7 C D^ix : First, the fixed point R|>, and second, the fixed point P. 
The latter is a saddle whose stable manifold (intersected with 3^ix) is two-dimensional; hence, there 
exists a one-parameter family of orbits that converge to P as r — > 00. The fixed point R^, however, 
is a sink, which means that there exists a two-parameter family of orbits that converge to R& as 
r — * 00. Accordingly, in the case A 1 *, there exists an open set of forever expanding cosmological 
models; the asymptotics of these models is governed by the approach to R|,; see Fig. 17. (The 
existence of recollapsing solutions is evident from the properties of the flow on X® x .) To find the 
past attractor of orbits in J^rx we follow the well-established principles that have already been 
successfully applied in the proof of Theorem 2: Analyzing Figs. 5, 13, and 15, it is immediate that 
the only possible a-limit set of typical orbits is the source Tt, cf. Table 6. This completes the proof 
of the lemma. □ 

Interpretation of Theorem 3. Cosmological models of Bianchi type IX with anisotropic matter 
satisfy the closed-universe-recollapse conjecture [2, 7, 21], i.e., these models recollapse. However, 
there is one restriction: The anisotropic matter model is required to be different from A 1 *, see (72) 
and Figs. 1, 10 and 11. (Recall that the energy conditions are satisfied for the A 1 * matter model.) 
In the case of anisotropic matter of that kind, type IX solutions need not necessarily recollapse; 
there exists a set of typical solutions (corresponding to an open set of initial data) that are forever 
expanding. We conclude that "anisotropic matter of the type A 1 * matters" (a lot). Let us thus 
discuss the anisotropy case A 1 ^ in detail: Toward the initial singularity, every typical solution with 
anisotropic matter of the type A 1 * is asymptotic to the Taub solution (30); see Fig. 17. However, 
there exist non-generic solutions that exhibit different kinds of behavior: We observe solutions with 
isotropic singularities, cf. (32), and solutions that approach (102) as t — > 0. Toward the future, 
recollapsing solutions and forever expanding solutions exist on an equal footing. Recollapsing 
solution exhibit a final singularity of the Taub type. Expanding solutions, on the other hand, are 
geodesically complete toward the future; the future asymptotics are characterized by an approach 
to the asymptotic metric (101). 



45 



Theorem 4. The qualitative behavior of typical orbits of the dynamical system (79) that represent 
the expanding phase of LRS Bianchi type IX models ( in the various anisotropy cases ) is the one 
sketched in Figs. 16 and Fig. 17. The past attractors for the various cases and the possible limit 
a-limit sets for non-generic solutions are given in Table 6. 

Remark. Theorem 3 states that, for every anisotropy case except A™, all models are recollapsing 
models. Theorem 4 discloses the possible past asymptotic behavior; the possible future asymptotic 
behavior of these recollapsing models arises from the past asymptotic behavior by applying the 
discrete symmetry (65). 

Proof. The proof of Theorem 4 follows the well-established principles that have already been 
successfully applied in the proof of Theorem 2. Lemma 2' restricts the possible a-limit sets to a 
subset of the boundary of the state space J^x- Since by the analysis of Sec. 10, the dynamics on 
these boundary subsets is fully understood, the proof of Theorem 4 essentially amounts to gluing 
together the phase portraits of Figs. 5, 13, and 15. 

The invariant structures on the boundary subsets (59) that are potential a-limit sets are fixed 
points, periodic orbits, and heteroclinic cycles/networks. Figs. 5, 13, and 15 show that periodic 
orbits do not occur, so that the a/w-limit sets of orbits in y^ must be fixed points or heteroclinic 
cycles/networks (if present). The investigation of the fixed points is straightforward. The procedure 
is completely analogous to the one in the proof of Theorem 2. The results of the local analysis of 
the fixed points is summarized in Table 6. Let us thus restrict ourselves to discussing the role of 
the heteroclinic cycles/networks. 

Heteroclinic cycles are present only in the cases D_, B_|_, C_|_, and D_|_, see Fig. 16. In the cases 
B_|_ and C_|_, inspection of Figs. 5 and 15 reveals that there does not exist any fixed point on 
iSj U X\ (cf. Lemma 2') that acts as a source; we simply use the color-coding of these figures. This 
leaves only one structure on <Sj U X \ that can be the a-limit set of typical orbits: the heteroclinic 
cycle depicted in Fig. 16(c). Case D+ is analogous; however, instead of a heteroclinic cycle there 
is a heteroclinic network, see Fig. 16(f). It is this network (or a subset thereof) that it the past 
attractor. Finally, consider the case D_ . In analogy with the arguments presented in the proof of 
Theorem 2 we find that the heteroclinic cycle see Fig. 16(a), is not a possible a-limit set of 
orbits of 3^rx- (The same is true for the cycles involving an orbit connecting Tj with Qj on X\.) 
This concludes the proof of the theorem. □ 

Interpretation of Theorem 4- Anisotropic matter of the types D_, C_, B_, A_\A™, and A_|_ "does 
not matter" . The asymptotic behavior toward the initial singularity of typical Bianchi type IX 
solutions (associated with anisotropic matter of one of these types) is the same as for typical perfect 
fluid solutions, see Sec. 5. We observe that the typical solutions behave like the Taub solution (30) 
as t — > 0. However, there exist non-generic solutions that exhibit different kinds of asymptotic 
behavior: There are solutions with isotropic singularities, cf. (32), solutions that approach (102), 
and solutions that approach (100) or (101) as t — > 0; see Tabic 6 for details. Toward the future 
we observe recollapse; the asymptotics toward the final singularity is analogous to the asymptotics 
toward the initial singularity. 

In contrast, anisotropic matter of the types B+, C+, and D+ "matters". Gcncrically, the approach 
of the associated LRS type IX models toward the (initial and final) singularity is oscillatory: The 
solutions oscillate between the Taub family (30) and the non-flat LRS family (31). This kind of 
asymptotic behavior differs considerably from that of (LRS) vacuum and perfect fluid solutions 
(as discussed in Sec. 5). 
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(d) A+ G8j| < p < 1) (e)B+,C+ (f) D+ 



Figure 16: Phase portraits of Bianchi type IX orbits in the various anisotropy cases — the special 
case A 1 * is depicted in Fig. 17. Bold lines are typical (in fact generic) orbits; dashed lines are orbits 
that lie on the boundary of the state space and connect to form heteroclinic cycles/networks. In 
the cases B_|_ and C_|_, the past attractor is the heteroclinic cycle; in case D_|_, the past attractor 
is the heteroclinic network (or a subset thereof). In the remaining cases, the past attractor is the 
fixed point Tj,. There exist non-generic orbits that display different asymptotic behavior toward 
the past, see Table 6. The phase portraits in the cases Ao-|- are the same as for A_. 

12 Examples of matter models 

In this section we present in detail two examples of matter models to which one can apply the 
main results of this paper: Ensembles of collisionlcss masslcss particles, described by the Vlasov 
equation, and magnetic fields. (For a discussion of a third matter model, elastic matter, we refer 
to [4].) Collisionlcss matter with masslcss particles satisfies our assumptions and the results of 
this paper apply directly. Magnetic fields violate Assumption 5, but we show below that our 
analysis extends straightforwardly to this matter model as well. In fact we believe that, up to 
minor changes, our method should work for an even larger class of matter models than the one 
considered in this paper. For a general introduction to collisionless matter and the Vlasov-Einstein 
system we refer to [1]; the Bianchi type I case is discussed in detail in [14]. For a generalization of 
the results presented here to Vlasov matter with massive particles we refer to [8] . 
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Figure 17: Phase portrait of typical LRS Bianchi type IX orbits in the special case A 1 — defined 
by (72). This is the only anisotropy case where expanding (i.e., non-recollapsing) solutions exist; 
these solutions satisfy Hd > for all times. The future attractor of these solutions is the fixed 
point R[> . Note that this type of behavior is typical. A second class of typical orbits corresponds 
to recollapsing models. Finally, there exist non-generic solutions whose w-limit set is the point P; 
these orbits are not depicted; instead we refer to Table 6. 



12.1 Collisionless matter 

Consider an ensemble of particles with mass m in the spacctimc (M, 4 g). This ensemble of particles 
is represented by a distribution function ('phase space density') / > 0, which is defined on the 
mass shell, i.e., on the subset of the tangent bundle given by 4 g(i>, v) = —m 2 , where v denotes the 
(future-directed) four momentum. If (£, x 1 ) is a system of coordinates on M such that dt is timelike 
and d x i is spacclike, then the spatial coordinates v l of the four-momentum are coordinates on the 
mass shell, and we can regard / as a function / = f(t, x 1 1 v J ), i, j = 1,2,3. The energy-momentum 
tensor is defined as 



The equation satisfied by the function / depends on the kind of interaction between the particles. 
If the interaction is through binary collisions, the resulting equation is the Boltzmann equation [12]. 
Here we assume that the particles are freely falling, i.e., they interact only through gravity. By the 
equivalence principle, the particles trajectories must coincide with the geodesies of the spacetimc. 
Accordingly, / has to be constant along the particles trajectories and thus solves the Vlasov 
equation 




(84) 
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Fixed point 



Dimension of unstable [stable] manifold intersected with 3^x 



A_ (not A?) 



A + (/3</3„) A + (/3 9 </3) 



B 



+ 



'+ 



R b 

F 

c s 

hct. cycle 
P 



3 

[3] 

2 
1 

[2] 



3 

2 



1 
3 



1 
3 



1 
3 



Table 6: Possible a-limit sets of LRS Bianchi type IX models and w-limit sets of forever expanding 
LRS Bianchi type IX models. Each orbit in J^x corresponds to (the expanding phase of) a Bianchi 
type IX solution (and conversely); hence the intersection of the unstable [stable] manifold of a fixed 
point with yields the set of type IX solutions converging to the solution represented by that 
fixed point as t — > [t — > oo] . If this set is three-dimensional, the fixed point is a source [sink] ; if 
the dimension is two, there exists a one-parameter set of orbits converging to the fixed point; if 
the dimension is one, there is only one orbit with that property. Center manifold analysis shows 
that the cases A -|_ and A _ behave like A_|_ (0 < (3 < f3$) and A_ (not A™), respectively. The 
exact solutions represented by the fixed points are given in Sec. 5 and in Appendix A. 



(The characteristics of the Vlasov equation — which correspond to the particles trajectories and 
along which / is constant — coincide with the lift of the geodesic flow on the tangent bundle, in 
accordance with the geometric interpretation of the free fall motion of the particles.) 

Consider now a spacetime of Bianchi class A with LRS Bianchi symmetry and define a time 
independent orthogonal frame such that the metric takes the form 

4 g = -dt 2 + .911 (t) LJl ® 0>1 + #22 (0(^2 <8> U>2 + ^3 ® w 3 ) ■ (86) 

As proved in [23], the general solution of the Vlasov equation (85) on a background spacetime with 
the metric (86) can be expressed as 

f = fo(vi,v 2 2 +vi), (87) 

where fo : R x K + — > M + is an arbitrary, sufficiently smooth function. The function fo can be 
interpreted as the 'initial data' for / at some time t = to- 

However, for the pair (86) and (87) to be a candidate for a solution of the Einstein- Vlasov system, 
the energy momentum tensor must be compatible with the structure of (86), i.e., it must be 
diagonal and satisfy T\ — T\. This can be achieved by restricting to distribution functions (87) 
that are invariant under the transformation v\ — > —v±. These distribution functions are called 
reflection symmetric, see [28]. We assume further that the support of fo docs not intersect any 
of the axes (split support assumption), which ensures that p > 0. It is now clear that, with the 
hypotheses on fo, Assumptions 1-3 are satisfied. 

Assumption 4 is only satisfied in the case of massless particles. For m = we have 

T f r I 11 2 , 22, 2 , 2^1/2 dvidV2dv 3 

P = T Q = J fo(g v 1 + g (v 2 + v 3 )) , (88a) 

n^i I r it 2 ( 11 2 , 22/ 2 2^-V2 dvidv 2 dv 3 

Pi = T l = J fo g 11 < V + 9 (v£ + « 3 )) ; (88b) 

note that P2 = P3 , because / is invariant under the exchange of «2 and v 3 . It follows that 

V = I (pi + Vi + Ps) = I P , 
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i.e., w = |. (For m > 0, the equation of state between the isotropic pressure and the energy 
density is non-linear.) The rescaled principal pressures are given by 

f , n „ v J fo [0- - 2s W + S ( V Z + "3 )] " V2 ^i 2 dvidv^ reQ , 

«Ji(s) = (1 - 2s) - — , (89a) 

/ f (v) [(1 - 2s)v* + + « 3 2 )] 1/2 d Vl d« 2 d« 3 

, s , x ,s / /o [(1 - 2g)t'! 2 + s{v 2 + « 3 2 )] " 1/2 vj d Vl dV2dv 3 
u{s) = w 2 {s) = w 3 (s) = s —r^—^ ;— — ; — — o „ 1/2 ; — ; — ; ■ (89b) 



/ f (v) [(1 - 2s)v? + s(v 2 + w 3 2 )] i/z d Vl dv 2 dv, 



3 



Therefore, using the terminology of Sees. 3 and 4, we obtain u(0) = and w(i) = i; hence 
Assumption 5 is satisfied with w_ = and v+ = 1. Moreover, a direct calculation shows that u(s) 
is monotonically increasing, see [28], so that Assumption 6 is satisfied as well (i.e., there exists a 
unique s such that u(s) = ^, whence w\ = w 2 = W3 = ^ at s). 

In conclusion, Vlasov matter with massless particles falls into the class of matter models considered 
in this paper. Since w = | and u_ = 0, we have (3 = 1, see (28), and therefore, Vlasov matter with 
massless particles is of type B_|_. The qualitative dynamics of solutions of the (Bianchi class A LRS) 
Einstein- Vlasov equations for the various Bianchi types is the one depicted in Figs. 2(g), 7(g), 16(e). 

Remark. If we consider ensembles of collisionless particles with positive mass, m > 0, w and the 
rescaled pressures Wi are functions of (s±, 52,33) and an additional scale variable, which can be 
taken to be z — m 2 /(m 2 + g 11 + 2g 22 ). In [8] we will show that the analysis of this paper carries 
over to this more general situation. 



12.2 Magnetic fields 



For an electromagnetic field represented by the antisymmetric electromagnetic field tensor F^, 
the energy-momentum tensor is given by 



TV = I 

4tt 1 



The equations for the field are the Maxwell equations 

, v 







In an LRS Bianchi class A spacetime with metric (86), consider specifically a purely magnetic field 
that is aligned along the axis perpendicular to the plane of rotational symmetry, i.e., 



F.j 



/0 0\ 



K 

\0 -K J 



the magnetic field is determined by K through B 1 = K(g 11 ) 1 ^ 2 g 22 , B 2 = B 3 = 0. The Maxwell 
equations imply that K is a constant; this implies that the energy density p, 



P 



8 71 



G? 22 ) 2 ^ 2 



is a function of the metric (which depends on the initial data for the magnetic field). Furthermore, 
T^ v is diagonal and 



T 



T . 



P- 



(90) 



Accordingly, Assumptions 1-3 are satisfied and diagonal models exist. It follows from (90) that 
-1 and u>2 = W3 = 1, so that w = ^. Thus Assumption 4 is satisfied, but Assumptions 5 



■50 



and 6 arc violated; in particular, there exists no isotropic state of the matter. Despite this fact, 
the analysis for magnetic fields cosmologies can be carried out as straightforwardly as the analysis 
for other matter fields by using the 'building blocks' discussed in this paper: 

Bianchi type I. The reduced dynamical system for LRS Bianchi type I aligned magnetic fields 
cosmologies is given by setting u — u>2 — if 3 = 1 and w — | in (42), i.e., 

S' + = (1 - S+)(2 - E+) , s' = -6s(l - 2s)S + . (91) 

Due to the decoupling of these equations, the original state space can be reduced to the one- 
dimensional state space E + £ [—1,1]. The equation for E + is identical to (46t,) with w = 1 and 
j3 = —2 (or, cquivalcntly, to (46j) with w — ^ and (3 = 4). Since E' + > 0, we conclude that the 
Taub point is the a-limit and the non-flat LRS point is the w-limit of all orbits. 

Bianchi type II. The reduced dynamical system for LRS Bianchi type II aligned magnetic fields 
cosmologies is derived from (49) by setting u = W2 = W3 = 1 and w = 3, i.e., 

E' + = {\Mf +0)(2 - E+) , s' = -6s(l - 2s)E+ , (92a) 

M[ = Mi [2(1 - 2E+) - \Ml - ti\ , (92b) 

where D, = 1 — Y? + — j^Mf. The equation for s decouples which implies that we may restrict 
ourselves to the reduced state space {-1 < £+ < 1, < Mi < [12(1 - £+)] 1/2 }• The system of 
equations on this state space coincides with the system (57) on the set S\, when we set w = | and 
[3 = —2. The phase portrait thus corresponds to the phase portrait of type D_ on S\,, see Fig. 6(a). 
(The system (53) on the set Sj with w = | and j3 = 4 is identical; see the phase portrait of type 
D_|_ on <Sj, Fig. 5(c).) Accordingly, the Taub point is the a-limit and the non-flat LRS point the 
kj-limit of every LRS Bianchi type II orbit. 

Bianchi type IX. The reduced dynamical system for LRS Bianchi type IX aligned magnetic fields 
cosmologies is obtained from (79) on by setting u = u>2 = w% = 1 and w = |, i.e., 

r = 2 r (H D (q - H D Y>+) - 3E+ sin 2 ■&) , (93a) 
& = -3E+ sin(2tf) , (93b) 
E' + = rsint?- 1 + (H D - E + ) 2 + H D Y, + (q - H D Y,+ ) + 2fl , (93c) 

where = 1 — Y? + — j r sin$, q = 2E^ + f2 and Hp = \/l — 2r cos d. The state space is J^x, 
see (78). (By using the formulation of Subsec. 10.2 we restrict ourselves to the dynamics of 
solutions in their expanding phase Hjj > 0; by applying the discrete symmetries (65) the entire 
dynamics is obtained.) The analysis of the system (93) is completely analogously to the analysis 
of Sees. 10 and 11. First, we analyze the flow induced on the four boundary subsets. On the side 
iSjj, the system induced by (93) corresponds to (69) with w = h and /? = 4, i.e., to the case D+ 
of Fig. 5(e). The dynamical system induced on the base B\x coincides with (71) for w = | and 
[3 = —2 (which is of type D_); see Fig. 13(a). The flow on the vacuum boundary is independent 
of the matter model. On the Bianchi type I boundary X\ the flow is equivalent to (91). Second, 
in analogy to the procedure of Sec. 11 we use the monotone function (83) and glue together the 
pieces. The final conclusion is that the fixed point Tj is the a-limit of every LRS Bianchi type IX 
orbit. Furthermore, every model is recollapsing and the w-limit set is the Taub point arising from 
Tjj by applying the discrete symmetry. 

Remark. The qualitative behavior of magnetic field cosmologies considered in this section is exactly 
the same as that of vacuum models, see Sec. 5. 

Remark. The conclusions are identical if we add an electric field parallel to the magnetic field. 
The Maxwell equations then show that E 1 = L{g 11 ) 1 / 2 g 22 and E 2 = E 3 = 0, where L = const. 
The energy density becomes 8irp = (g 22 ) 2 (K 2 + L 2 ), so that the energy-momentum tensor is a 
functional of the metric; accordingly, (90) and its consequences remain valid. 
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Remark (Erratum). The statement in [5] that magnetic fields are a matter model of type D_ in 
the context of Bianchi type I is erroneous. The content of this subsection corrects this claim: 
Magnetic fields do not directly fit into the classification of Table 2 and Fig. 1, but can be treated 
in an analogous manner. 

13 Discussion and open problems 

In this paper we have studied in detail spatially homogeneous locally rotationally symmetric so- 
lutions of the Einstein-matter equations. Rather than explicitly specifying a particular matter 
source, we have imposed a set of assumptions on the matter that characterize a large class of 
matter models including classical examples like elastic matter (associated with a large variety of 
equations of state), collisionless matter (Vlasov matter), and magnetic fields; in addition, the per- 
fect fluid matter model with a linear equation of state, which is the most commonly used matter 
model in cosmology and astrophysics, is naturally embedded in the class of matter models we 
consider. The main aim of our analysis has been to study the influence of the matter source on the 
dynamics of the associated (SH LRS) cosmological models. A natural focus has been to ask how 
robust the qualitative dynamics of perfect fluid solutions is under changes of the matter model. 
Does matter 'matter'? 

The main result of our analysis is that, indeed, anisotropies of the matter model 'matter'. Two 
facts arc of particular interest: (i) There are anisotropic matter models that satisfy the strong (and 
the dominant) energy condition such that the associated cosmological solutions of Bianchi type IX 
(i.e., closed cosmological models) do not recollapse but expand forever; this is in stark contrast to 
perfect fluid models and other anisotropic matter models; see Theorem 3 and its interpretation in 
Sec. 11. (ii) For a class of matter models including collisionless matter the asymptotic behavior 
of solutions toward the initial singularity is completely different from that of vacuum and perfect 
fluid solutions in that the approach to the singularity is not governed by the Taub solution but 
oscillatory; see Theorem 4 and its interpretation in Sec. 11. 

There are several interesting open problems. We have restricted our attention to three Bianchi 
types of class A: Bianchi type I, II and IX. However, as already indicated in Sees. 6 and 9, it 
is clear that the methods we have developed can be extended and/or modified to analyze other 
cosmological models of interest: On the one hand, there are the Kantowski-Sachs models, which 
are interesting because their spatial topology is S 1 x S 2 and thus closed. (Do models necessarily 
recollapse?) On the other hand, there are the Bianchi type III models (which are of class B) and 
the type VIII models, the latter being as fundamental as type IX models for our understanding 
of generic spacelike singularities. Note that in the analysis of type III and type VIII models an 
additional challenge might arise, since in each case the state space for the reduced dynamical 
system describing these models is probably unbounded. We leave the analysis of the Kantowski- 
Sachs models and the LRS Bianchi type III and type VIII models as interesting open problems. 

We expect that there exists a number of ways to extend the analysis and the results presented 
in this paper. Let us begin with Assumption 1. There exists an obvious generalization of this 
assumption, which is to require that the stress-energy tensor is a function not only of the spatial 
metric but of the second fundamental form fc.y as well. This new assumption would lead to a 
generalization of eq. (9) and, in the LRS case, give rise to rescaled principal pressures that depend 
not only on s but also on S + . 

Let us next turn our attention to Assumption 4. This assumption imposes a linear equation 
of state between the isotropic pressure and the energy density, which implies that p depends 
on n = 1 / Vdet g as indicated in (16). However, based on the results derived in this paper, 
it is relatively simple to treat nonlinear equations of state. Let us elaborate. Assume that p 
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is given by a more general function p = p(n, si, S2, S3) instead of (16). An interesting subcase 
is p(n, s\, S2, S3) = ip(n)ip(si, S2, S3). In the latter case we obtain that Wi = iUj(n, Si, S2, S3), 
i = 1,2,3, satisfy (17) where w = w(n) = (d\ogip/dn) — 1. In this case, the Einstein evolution 
equations, written in terms of the dynamical systems variables of Sec. 6, decouple into an equation 
for D and a reduced system of equations for the remaining variables, which is given by (37), 
supplemented by an additional equation for n, i.e., n' = — in. Assume we have fixed the Bianchi 
type and let X denote the state space for the corresponding reduced dynamical system. When we 
choose to compactify the variable n, i.e., when we replace n by N = n/(l + n), the state space of 
this dynamical system is X x (0, 1). If we assume an equation of state such that Wi(N, s±, S2, S3) 
possess well-defined limits as N — > and N — > 1, we can extend the dynamical system to the 
boundaries Xq = X x {0} and X\ = X x {1} of the state space. The dynamical system on each 
of these boundary subsets coincides with the system (37) that we have discussed so extensively 
in this paper. In particular, we may define a parameter (3q using the rescaled principal pressures 
Wi(0, Si, S2, S3) in (19) and a parameter /?i using Wi(l, si, S2, S3) and apply our classification of 
matter models for N = and N = 1 separately. Of course, in general, the matter model on 
the boundaries TV = and N = 1 will not be of the same type. Since the variable N is strictly 
monotone, the asymptotic dynamics of solutions of the dynamical system is associated with the 
limits N — > and N — > 1. Accordingly, asymptotically, the flow of the boundary subsets Xq and X\ 
(and thus the results of the present paper) constitute the key to an understanding of the dynamics 
of the more general problem with nonlinear equations of state. 

Finally, consider generalizations of Assumptions 5 and 6. Concerning the latter, it is possible 
to permit the existence of more than one isotropic state or the case that isotropic states do not 
exist at all. The latter is satisfied by magnetic fields, see Subsec. 12.2, which thus seem more 
interesting from the physical point of view. A natural way to generalize Assumption 5 is to require 
the existence of three constants c\, C2, C3, such that Wi(s\, S2, S3) = Ci when s,; = 0. A simple 
subcase is Wi(si, S2, S3) = q, i.e., the rescaled pressures are constants, which holds for magnetic 
fields, see Subsec. 12.2. If we generalize Assumption 5 as indicated, there would appear three 
parameters, instead of /3 alone, in the reduced dynamical system. In this case the number of 
possible dynamics and bifurcations would increase considerably, but no conceptual new difficulty 
would arise. Whether new phenomena occur under these relaxed assumptions is an open question. 



A Exact solutions 

In the course of our analysis of the system of equations (37) we have discovered a number of fixed 
points, see Table 7. In this appendix we derive the exact solutions that correspond to the these 
fixed points. 

Consider an arbitrary fixed point of Table 7 and let (s, £+, Mi, Hd) be the coordinates of this 
fixed point; the associated deceleration parameter is 

q = 2T? + + \{l + iw){l-T 1 \-^n\Ml) , 
see (38) and (39'). Expressed in i, eq. (36) reads 

D- 2 d t D = -[H D (l+q) + ^ + (l-H 2 D )}, (94) 
which can be solved to obtain 

D(t)=[H D (l + q) + j: + (l-H 2 D )]- 1 t- 1 . (95) 

where we have shifted the origin of time so that D diverges as t — > 0. Using (35b) and (39') this 
leads to 

pit) = i ± — 12 1 1J „ r 2 . 96 

[H D (l + q ) + Z + (l-H 2 D )} 2 
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Fixed 










Solution 


Bianchi 


point 


s 




iw l 


type 


T b 





-1 





1 


(30) 




Tj 


1 

2 


1 

1 





1 
1 


(30) 


1 




U 


1 
1 


U 


i 
1 


(31) 




Qb 


1 

2 


1 





1 


(31) 




F 


S 








1 


(32) 




Rb 





-fi 





1 


(101) 




Rj 


i 

2 


P 
2 





1 


(100) 




c b 





l+3«j 


36(l-w)[3/3 2 (l— tu)+8/3+(l+3i«)] 


1 


(103) 


II 


8+3/3(1— tu) 


[8+3/3(1-™)]* 




1 


l+3u> 


36(l-tu)[3/3 2 (l-M>)-16/3+4(l+3uj)] 


1 


(102) 


II 


2 


16-3/3(1— uj) 


[16-3,3(1— uj)] 2 


p 





l+3tu 





2+3/3(1— w) 


(104) 


III/KS 


^/(l-3u;) 2 +6/3(l— tu) 


^/(l-3iu) 2 +6/3(l— uj) 



Tabic 7: Tabic of the fixed points of the reduced dynamical system. Each of the points R;,, Rj, Ct,, 
Cj, and P exists only under certain restrictions on the anisotropy parameter j3. 



The principal pressures pi(t), p 2 (t) are then given by 

Pi(t)=Wip(t) = (3w-2u(a))p, p 2 {t) = w 2 p = u(s)p , (97) 

see (27), where the constant u(s) is known (since s G {0, s, h}), see Sec. 4. 
To obtain the metric we note that 

d mi = -2 9ll D(t) (2E+ - H D ) , d t g 22 = 2g 22 D(t) (E + + H D ) , (98) 

which follows from (12b) in combination with (22) and (35a). Inserting (95) into (98) and solving 
we obtain 

9ll = at 2 ^, g 22 = bt 2 ^, (99a) 
where a, b are positive constants and 

H D - 2£ + + 
71 ff fl (l + ?) + E + (l-fl|,) ' 72 F D (l + e?) + S + (l-if£,) ' lJyDj 

We remark that the constants a, b are in general not arbitrary: Imposing the Hamiltonian constraint 
may give rise to a restriction on their values. 

Applying the above algorithm to the fixed points of Table 7 yields the exact solutions represented 
by these points. For T|,, Tj we obtain the Taub solution (30); for Qi,, Qjj we obtain the non-flat 
LRS Kasner solution (31); for F we get the Robertson- Walker metric (32). 

The fixed points Rj and Rt, are of Bianchi type I, i.e., hi = and h 2 = h% = 0. 

Fixed point Rj : 

71 = 8 ALl3 72 = W + fi (100 a) 

71 3[/3 2 (l-w)+4(l + w)] ' 72 3[/3 2 (l-w)+4(l + u;)] ' ( ' 

16(4 - 3 2 )t~ 2 

P = 3[^(1- W )+4(1 + W )] 2 ' Pl = X 2 + & ^ > P 2 = * - ® + ^ ■ ( 100b ) 
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Figure 18: The values of 71 and 72 for the exact solution (101) corresponding to R\,. Each curve 
corresponds to a given value of w, where the topmost curves correspond to w = — | and the 
lowermost curves to w = 1. 



Fixed point Rt, : 

2(1 + 2/3) 2(1-/3) 
ri 3[/3 2 (l-w) + l + u>] ' 72 3[/3 2 (l-™) + l + w] ' 1 ' 



P = 7 



4(l-/3 2 )i" 2 



3[/3 2 (l-u.) + l + ' 



Pl = Kl-/?)+/3]p, p 2 = i[ W (2 + / 3)-/3]p. (101b) 



Let us analyze the solution (101) associated with R^. The values of 71 and 72 vary considerably 
depending on w and on the anisotropy parameter (3 (where the existence of Rt, requires — 1 < (3 < 1). 
In Fig. 18 we depict 71 and 72 as functions of (3 where w takes a fixed value. The topmost curves 
corresponds to the value w = — ^; the lowermost curves correspond to w = 1. An interesting 
observation is the possible occurrence of partial (directional) accelerated expansion for the metric 
given by (101). It is straightforward to see that 71 > 1 is possible for some range of /3 provided 
that w < ^. (The range is (3 € (0, 1) for w = — A and decreases to [3 £ (1 — e, 1) when w approaches 
w = i.) The maximal rate of acceleration is obtained by maximizing -fi over the admissible domain 
of (3 and w; this yields the maximal value of 71 = k(1 + \/3) ~ 1.366, which can only be attained 
in the limit w -> -1/3 and p -»• \(-\ + w 0.366. 

Analogously, a straightforward calculation reveals that 72 > 1 is possible, which means that lengths 
in the plane of local rotational symmetry expand at an accelerated rate, see Fig. 18. The condition 
for 72 > 1 is w < i(l — -v/3) and [3 £ (/?-, (3+), cf. (72b), i.e., 72 > 1 in the anisotropy case A 1 *. The 
maximal rate of acceleration is obtained by maximizing 72 over the admissible domain (72); this 
yields the maximal value of 72 ~ 1.112. This value can only be attained for w close to —1/3; for 
larger values of w the maximal acceleration is lower (e.g., the maximum of 72 is 1.067 for w = —0.3 
and 1.007 for w = —0.25). While expansion of areas might be accelerating, the standard volume 
expansion is decelerating; we have y/det g oc i 3 /( 1+| z) hence the length scale behaves like t 1 /( 1+< ?) . 

The fixed points Cj and Cj> are of Bianchi type II, i.e., hi = 1 and n 2 = h 3 = 0. These fixed points 
correspond to solutions that naturally generalize the Collins-Stewart solution (33). 

Fixed point Cj : 

(4-/3)(l- W ) _ 2/3(l + w )+6-/3 

71 8(1 + w) - /3(1 - w) ' 72 8k;(1 + ff) + 8 - /3 ' [ ' 



1 -j£rr^ W)) w ■ Pi = \[w{2 + f3)-P] P , p 2 = \[w(4-(3) + f3}p. (102b) 

[8(1 + w) — p(l — w)Y 
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Furthermore, the constants a, b of (99) are related by 



y/E_ y/(l- w)[3^(l -w)- 16/3 + 4(1 + 3w)} 

2b ~ 8 - /3 + 8^(1 + /3) ' 1 > 

where the function under the square root is positive for ft < /3j (which is the prerequisite for the 
existence of Cj). 

Fixed point Ct,: 

(2 + /?)(! -a;) _ 3 + /3-2/3(l+ W ) 

71 4(1 + w)+ /3(1 -w) ' 72 4w(l - 2/3) + (3 + 4 ' 1 ' 



P= 8 f An Z\ 2 1 [ L W)) V > Pi = [w(l-ft)+ft]p, p 2 = ±[w(2 + ft)-ft]p. (103b) 
[4(1 + u>) + /3(1 — u>)J z 



Furthermore, a and b are related by 



^ _ y/(l - H[3/? 2 (l - w) + 8ft + (1 + 3w)] 

2b ~ 4 + /3 + 4^(1-2/3) ' 1 J 

where the function under the square root is positive for (3 > ft, (which is the prerequisite for the 
existence of Ct). Note that passing from Rjj to R\, and from Cj to C\, corresponds to replacing (3 
by -2/3. 

Fixed point P: 

(3(1 -w)-2w . > 

7i = r , 72 = 1 , 104a 

1 + w + /3(1 — w) 

P= u 2 l {1 l^\~\2 r ^ Pi = [w(l - (3) + ft]p , P2= i[ w (2 + ft)-p]p. (104b) 
[1 + w + /3(1 — w)^ A 



The fixed point P is a fixed point on the Kantowski-Sachs subset Six of Xix] its existence is 
restricted to the case A™, which is defined through (72). The point P represents a solution of 
Kantowski-Sachs type, cf. (116); see Sec. 9.2 and the remark in Subsec. 10.1 for details. The 
constant b of (99) is given by 

b = _(l- w)m i- W ) + 2ft + (l + 3 W )] ; 
(1 + w + J3(l - w)) 

since the numerator is negative, cf. (73), b is positive as required by (99). 

The value of 71, see (104a), depends on w and the anisotropic parameter (3 (which are subject 
to (72), since the anisotropy case is A 1 ^). We find that 71 < 1 irrespective of w and /3; the possible 
values of 71 are depicted in Fig. 19. Accordingly, the S 1 component of the Kantowski-Sachs 
metric (116) expands at a lesser rate than the S 2 component. 

Remark. There exists a different range of the parameters w and ft (different from A 1 ^) such that P 
is located on the (type III) boundary of the Bianchi type VIII state space. In this case, the fixed 
point P represents a type III solution and the constant 6 of (99) is given by 

_ (1 - w) [3/3 2 (l - w) + 2/3 + (1 + 3w)] 
(1 + w + /3(1 - w)) 2 

which is then positive. 
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Figure 19: The value of 71 for the exact solution (104) corresponding to P. The anisotropy case is 
the Bianchi type IX special case A™, which is defined by (72); see Fig. 10. For a given value of w, 
the value of 71 depends on j3; if j3 =/?-, then 71 takes the value on the solid black line; if = (3+, 
then 71 takes the value on the dashed black line; if /?_ < j3 < 0+, then the value of 71 is in the 
gray region. The axes of the diagram are 71 = and w = —0.23. 



B Local rotational symmetry 

In this appendix we discuss in detail the concept of locally rotational symmetry (LRS). We analyze 
the LRS Bianchi class A models, as well as the Kantowski-Sachs models and the Bianchi type III 
models. 

Spatially homogeneous cosmological models that arc locally rotationally symmetric (LRS) are 
spacetimes that admit a four-dimensional isometry group whose orbits are three-dimensional space- 
like hypcrsurfaccs. We distinguish two cases: For the LRS Bianchi models the four-dimensional 
isometry group admits a three-dimensional subgroup Q that acts simply transitively on the three- 
dimensional orbits. For the Kantowski-Sachs models there does not exist any three-dimensional 
subgroup with that property. 



B.l LRS Bianchi class A models 

For a Bianchi model of class A there exists an adapted left-invariant frame £j on Q such that 

for any triple (ijk), see Sec. 2. (Throughout this appendix we employ the Einstein summation con- 
vention.) The structure constants are thus represented by a triple (ni, h%, n 3 ) with n, G {0, —1, +1} 
Vi, see Table 1 of Sec. 2. The spatial metric of a Bianchi model reads gij Cj 1 ® uP, where the com- 
ponents gij are (spatial) constants (i.e., functions of the time variable alone) and Co 1 is dual to 

For an LRS Bianchi model the isometry group of the spatial metric is four-dimensional. Therefore, 
there exists a one-dimensional isotropy group corresponding to a pointwise (left-invariant) symme- 
try of the metric. Visualizing this symmetry as a symmetry of the 'unit sphere' {X 1 1 gijX l X : > = 1} 
we see that merely one continuous isometry comes into question: Axial symmetry. An axially sym- 
metric 'unit sphere' is characterized by two degrees of freedom (represented by the axes of the 
ellipsoid); we conclude that the assumption of a four-dimensional isometry group will reduce the 
number of degrees of freedom of the spatial metric to two. 
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In this section we analyze under which conditions LRS Bianchi models exist — note that there are 
Bianchi types that are incompatible with the assumption of LRS symmetry — and we discuss their 
main properties. We begin with the simplest (and most illustrative) example. 

The simplest of Bianchi models are those of type I. Since 

[11,61 = 0, [6,6] = o, [6,6] = o, (105) 

the symmetry-adapted frame is in fact a coordinate frame; we denote the coordinates by x l . In the 
LRS case there exists a one-parameter isotropy group corresponding to a rotational symmetry 
w.l.o.g. we assume that the axis is 6 = d x \\ we thus have 6 i— > 6 and £/ i— ► 0^6f where 
O G SO (2); Latin capitals /, J, . . . run over 2 and 3 (and the Einstein summation convention is 
understood). The generator of this rotation in the 2-3-plane is 

3 <9 9 9 o ; 

Accordingly, the four-dimensional Lie group that represents the isometry group of an LRS type I 
Bianchi model is generated by the Lie algebra 

[6,6] =o, [6,6] =o, [6,6] = o, [f),£i] = o, [r?,6] = -6, [*?,6]=6- (106) 

Any metric that is invariant under the action of the Lie group generated by (106) must satisfy 
giK oc Sik, i.e., 

.922 = .933 ■ 

Remark. There exist several four-dimensional Lie algebras that contain (105) as a subalgebra, see, 
e.g., [27]; note, however, that (106) is the only Lie algebra that generates a Lie group that can act 
as an isometry group on a Bianchi type I spacetime. 

Let erj denote the standard permutation symbol (for the indices I, J = 2,3). For an arbitrary 
Bianchi type (of class A) we may write the commutator relations as 

[6, 6] = e LJ nx 6 , [6,6f] = B'Jj , where B = (j^ ~™ 3 ) . (107) 

Consider a rotation in 6, 6 with axis £i, i.e., let O G SO(2), and define = £i and 6 = K j£k- 
The tensor ejj is invariant under this transformation, because K jO L j€kl = £u- Hence the 
commutators (107) result in 

[i'n i'j] = eu ni & , [&,&] = (O^B O) 1 K & . (107') 
It is straightforward to see that B is invariant, i.e., 0~ 1 B O = B for all O G SO(2), if and only if 

n 2 = n 3 . (108a) 

The assumption 

.922 = .933 (108b) 
then yields an invariant metric, i.e., a LRS Bianchi class A model. 

Therefore, a spacetime of Bianchi class A admits a four-dimensional isometry group, if and only 
two structure constants are equal (which we assume to be hi — ns) and if the corresponding metric 
components are equal as well (i.e., gn = 333). The four-dimensional isometry group is generated 
by the Lie algebra 

[6,6] ="36, [6,6] =^16, [6,6] = ^26, (109a) 
W,6] = o, [j? ! 6] = -e"6, [*/,6] = e6, (109b) 
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where fi2 = fi3 and e = 1. Let us discuss the Bianchi class A models type by type. 

In the Bianchi type I case there exist three different but equivalent representations of LRS solutions; 
we make the choice (722 = 533- In the Bianchi type II case, for a fixed triple (711,712,713) = (1, 0, 0), 
there exists a unique class of LRS models, which is characterized by c/22 = 533 (since n 2 = 713 = 0). 

The Bianchi type VIo case is exceptional: There do not exist any LRS type VIo models. The reason 
is that the structure constants are pairwise different, which entails that the rotations O £ SO(2) 
do not act as isometries; note that this is irrespective of the assumption 522 = 533- Alternatively, 
we note that (109) does not define a Lie algebra in the type VIo case; setting (711,712,713) = 
(0, 1, -1) in (109) we find [fj, [|i,| 2 ]] + [| x> [| 2 , fj]] + [| 2 , [fj, |i]] = -e[r), | 3 ] + e[|i,| 3 ] = -2e| 2 + 0. 
However, while Bianchi type VIo is incompatible with a four-dimensional isomctry group, instead 
of a continuous symmetry there exists a discrete symmetry. The reflection 

o=(! J) (HO) 

leaves the matrix B in (107) invariant, i.e., 0~ 1 BO = B. The tensor eu is not invariant, but 
77,1 = 0; therefore, (107) is invariant under (110). The transformation (110) is a discrete isometry 
if 

.922 = 333 ■ 

Therefore, among Bianchi type VIo models there do not exist any LRS models; however, there exist 
models that admit a discrete isometry in addition to the three-dimensional group of isometries. 
Remark. In general, if 77,2 = —77,3, we find that the transformations O G SO(l, 1) define a continuous 
isometry for a metric with (722 = — 533- Hence, if the spatial metric were not a Riemannian but 
a Lorentzian metric, the models of Bianchi type VIo would admit a natural subclass of models 
whose isometry group is four-dimensional; one might call these models "locally boost symmetric" . 
For the Lie algebra of this isometry group we would have [77, £2] = ££3 and [77, £3] = ££2 in (109b). 

LRS models of Bianchi type VIIo admit an isometry group generated by 

[€1,6] = Is, [la,ls] = 0, [13,1x1= la, W.li] = 0, bUa] = -| 3 , ft, Is] = la- (HI) 
Define 

li = l(li+*7), l 2 =la, l 3 = ls, = 5 Hi +*/■)• 

We find that 

[li,la] = 0, [la,l 3 ]=0, [i' 3 ,i'i]=0, ftMi] = 0, [r?',| 2 ] = -| 3 , [#,&]=&, (111') 

hence the Lie algebra (111) coincides with the Lie algebra generating the isometry group of LRS 
type I models. The action of the subgroup generated by {^1,^2,^3} must be simply transitive, 
see [11, Appendix B]. Therefore, LRS Bianchi type VIIo models are in fact of type I. 

Finally, there exist LRS models of Bianchi type VIII and IX; in the latter case there exist three 
different but equivalent representations of LRS solutions since hi — hi — 77,3; we make the choice 

.922 = #33- 

Under the assumption of local rotational symmetry, where we make the choice (108), i.e., 77 2 = 773, 
cf. Table 1, and 1722 = 533, we find that the spatial Ricci curvature (11a) becomes 

R 1 1 = \ h\ m\ , R 2 2 = 77,177-2 7771777-2 — | 77^ 777^ , (112a) 

where 

777l = — , 7772 = Vg U ,9 U = 777-2 , 9^ = ^1^2 ■ (112b) 

The (spatial) curvature scalar (lib) reads 

i? = 77-1 777i + 2 77-1772 77717772 . (H3) 
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B.2 Kantowski-Sachs models 



For Kantowski-Sachs models, by assumption, the orbits of each three-dimensional subgroup G3 of 
the isometry group Q4 are merely two-dimensional. (Note that a four-parameter Lie group necessar- 
ily admits a three-dimensional subgroup, see, e.g., [11, Appendix A].) A two-dimensional manifold 
admitting a three-parameter isometry group is necessarily a manifold of constant curvature; this 
entails that the generators £j of the group G3 form the Lie algebra 

[£1,6] =6, [h,ta] = kii, (H4) 

where k S {0, — 1, +1} denotes the normalized curvature of the constant curvature orbits. (These 
Lie algebras are of Bianchi type VIIo, VIII, and IX, respectively.) However, in the cases k = 
and k = — 1, the existence of a three-dimensional subgroup G3 according to (114) implies that 
the Lie group G4 possesses additional three-dimensional subgroups, whose orbits are in fact three- 
dimensional, see [11, Appendix B] and [19]; accordingly, k = and k = — 1 arc LRS Bianchi 
models. 

On the other hand, in the case k = +1, if (114) is a subalgebra of a four-dimensional Lie algebra, 
then that Lie algebra is necessarily represented by 

[li,&]=&, li, [ls,|i]=6», (H5a) 

[V,ii]=0, K>6] = 0, [j?,|3]=0, (115b) 

see [11, Appendix B]. The Lie group Q± generated by (115) possesses a unique three-dimensional 
subgroup, the group G3 that is generated by (115a), cf. [11]. Therefore, it is possible to assume 
that C/4 acts as an isometry group of the spacetime in such a way that the orbits of G3 are two- 
dimensional spaces of positive constant curvature, i.e., 2-spheres. As a consequence, the metric 
reads 

4 g = -dt 2 +g n (t)dr 2 +g 22 (t) 2 g [s2] , (116) 

where 2 g[$2] is the standard metric on the 2-sphcrc. These are the Kantowski-Sachs models. 
Assuming that the coordinate r ranges in S 1 , the spatial topology is S 1 x S 2 and thus compact. 

The Ricci curvature Rij of the spatial part of the Kantowski-Sachs metric (116) has a rather simple 
structure. As suggested by (116) let x 1 = r and x 1 , 7 = 2,3, be coordinates on the 2-sphere. Then 
Ru = 0, Rik = 2 Sik (where I,K = 2,3), and the remaining components of the Ricci tensor 
vanish. Accordingly, 

R\ = , R 2 2 = g 22 = mim 2 , (117) 

where we resort to the quantities mi, m 2 of (112b). The component R 3 3 is identical to R 2 2 ; the 
remaining components are zero. The curvature scalar R is simply R = 2m\m 2 . 

Remark. The Ricci tensor (117) of Kantowski-Sachs models is obtained from the Ricci tensor (112) 
of LRS Bianchi class A models by formally setting n\ = and hxn 2 = 1. 

B.3 Lie contractions 

Contractions of a Lie algebra are obtained by considering sequences of basis transformations whose 
limit is singular 'in a controlled way', i.e., in the limit one observes convergence of the structure 
constants [31]. In [18] the notion of Lie algebra contractions is applied to the Bianchi models: 
For class A models there exists a hierarchy of Lie algebra contractions, which corresponds to 
successively setting the structure constants hi, h 2 , 77.3 to zero. In this way, the type IX algebra 
generates the algebras associated with the 'lower' Bianchi types VIIo, II, and I, and the type VIII 
algebra generates the algebras of the 'lower' types VIIo, VIo, II, and I. This hierarchy of the Bianchi 
types is of fundamental importance in the analysis of the dynamics of the 'higher' Bianchi types 
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- the asymptotic dynamics of 'higher' Bianchi types are directly related to the dynamics of the 
'lower' types, see [18] or, e.g., [33]. 

Consider the Lie algebra (109) of the four-dimensional isometry group of LRS type IX models, i.e., 
h\ = 1, h 2 = n-3 = 1 and e = 1. There are three possible contractions: We may set hi = 0, or 
h 2 = U3 = 0, or e = 0. Setting hi = (while leaving the remaining constants unchanged) we obtain 
the Lie algebra of LRS type VIIo models; this algebra coincides with the LRS type I algebra, see 
Sec. B.l. Setting h 2 = h% = we obtain the Lie algebra of LRS type II models. Finally, setting 
e = we obtain the Lie algebra (115) representing the isometries of Kantowski-Sachs models. For 
a detailed discussion of the contraction of (109) see [31]. 

Remark. The type IX Lie algebra (109), where hi = 1, hi = ^3 = 1 and e = 1, is isomorphic to the 
Lie algebra (115). In other words, the Lie algebra contraction obtained by setting e to zero does 
not yield anything new on the level of the Lie algebras (or the associated Lie groups). However, 
when (109) is regarded as an isometry group (on an LRS type IX space), the Lie algebra contraction 
process involves a contraction of the action of the group. In the singular limit that corresponds to 
the Lie algebra contraction one obtains a different action of (109); this action corresponds to the 
action of (115) on Kantowski-Sachs spaces. 

B.4 LRS Bianchi type III 

The Lie algebra of the four-dimensional isometry group of LRS type VIII models is given by (109) 
with Hi = — 1, h 2 = hs = 1 and e = 1. There are three possible contractions: We may set hi = 0, 
or hi — h% — 0, or e = 0. Setting hi — (while leaving the remaining constants unchanged) 
we obtain the Lie algebra of LRS type VIIo models; this algebra coincides with the LRS type I 
algebra, see Sec. B.l. Setting h 2 = ^3 = we obtain the Lie algebra of LRS type II models (in a 
different representation where the sign of hi is negative instead of positive). Finally, setting e = 
we obtain the Lie algebra (115) 

[li,6]=| 3 , [6,1a] = -&, [| 3 ,a]=6, (H8a) 
W,li]=0, 0,61 = 0, [*),&]= 0. (118b) 

The remarks of Sec. B.3 apply analogously. In [11, Appendix B] it is shown that this Lie algebra 
contains not only the subalgebra (118a) but also a subalgebra of the type 

= 0, = [&£l = o. (us) 

This is a Lie algebra of Bianchi type III, see, e.g., [33]. 

The action of the Lie group generated by (118) resembles the action of (115) on the Kantowski- 
Sachs models. Consider the metric 

4 g = -dt 2 + gn(t) dr 2 + g 22 (t) 2 g [H 2 } , (120) 

where 2 g[H 2 ] is the metric of constant negative curvature on hyperbolic space. This space is an 
orbit of the subgroup generated by (118a), i.e., (118a) acts multiply transitively on 2 g[H 2 ]- 

In analogy to (117), the Ricci curvature Rij of the spatial part of the LRS type III metric (120) 
has a rather simple structure. We obtain 

R\ = , R 2 2 = -g 22 = -mum , (121) 

where we resort to the quantities mi, m 2 of (112b). The component i? 3 3 is identical to R 2 2 ; the 
remaining components are zero. The curvature scalar R is simply R = 2mim 2 . 

Remark. The Ricci tensor (121) of LRS Bianchi type III models is obtained from the Ricci ten- 
sor (112) of LRS Bianchi class A models by formally setting h\—Q and hih 2 = — 1. 
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C Dynamical systems toolkit 



In this appendix we briefly present some concepts from dynamical systems theory that are used in 
this paper. For a thorough development of these concepts we refer to the book [26]; for a discussion 
of applications in cosmology we refer to [9, 33]. 

Let / : R" — * R" be a C 1 vector field and consider the autonomous dynamical system 

x = f(x) (122) 

for the function x = x(t) £ R". We assume that for every initial data point xq, the system (122) 
admits a unique global solution x £ C X (R), x(0) = xq. A point x* £ R™ is said to be an a-limit 
point of the solution x(t) if there exists a sequence of times t n — ► — oo such that x(t n ) — > as 
n — ► oo. The concept of w-limit is defined analogously by considering a sequence of times t n — > +oo. 
The a-limit set of a solution is the set of all its a-limit points. By the autonomy property, an orbit 
7 in the state space is associated with a one-parameter set of solutions; if the orbit reduces to a 
point, the associated solution of (122) is a stationary solution. Since solutions with the same orbit 
7 have the same limit sets, we may speak of a- and w-limit sets of orbits and we shall denote these 
sets as 01(7) and 0^(7). The following lemma collects some properties of a(j) that are used in the 
paper; analogous properties hold for 0^(7). 

Lemma 3. Let 7 be an orbit of the dynamical system (122). The following holds: 

(i) 01(7) is closed. 

(ii) If the exists a compact set E C R™ and to £ (—00, +00] such that j(t) £ E, for all t < to, 
then 01(7) is a non-empty, connected, compact subset of E. 

(Hi) If P £ ct("f), then the whole orbit through the point P belongs to 0(7). In particular, the 
a-limit set is flow invariant. 

(iv) 1/01(7) consists of a point P only, then P is a fixed point. 

(v) A source P is the a-limit of all orbits in a neighborhood of P. If P £ a("f), then a(j) = {P}. 

Let us recall the definition of a hcteroclinic cycle. A connecting orbit 7 is an orbit for which both the 
a- and w-limit sets consist of a point (which is then necessarily a fixed point). If a(j) = ^(7) = {P}, 
the set {7} U {P} is called a homoclinic cycle. If 71, 72 are two connecting orbits such that 
{Pi} = a; (71) = 0(72) and P2 = ^(72) = c*(7i), then the set {71} U {72} U {Pi} U {P2} is called a 
hetcroclinic cycle. Likewise one can define heteroclinic cycles consisting of an arbitrary number of 
connecting orbits. 

In the paper we often make use of the monotonicity principle. It can be stated as follows (for an 
extension see [13]): 

Theorem 5. Let S C R™ be an invariant subset of the dynamical system (122) and Z : S — * R be 
strictly monotonically decreasing along the flow; set a = mi{Z(y), y £ S} and b = sup{Z(y), y £ S}. 
Let 7 be any orbit in S. Then 

aM C{ S £dS : lim Z(y) £ a} , wM C {s £ dS : lim Z(y) £ b} . 

y—>s y—>s 

For a dynamical system over a two-dimensional state space, there exist further results that help 
studying the asymptotic behavior of orbits without the need for a monotone function (which is 
often hard to find). The most useful result is the Poincare-Bcndixson theorem, which states the 
following: 
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Theorem 6. Consider the dynamical system (122) in the plane, i.e., n — 2, and assume that 
there are at most a finite number of equilibrium points. Then, for any orbit, the a-limit set (u- 
limit set) can only be one of the following: A fixed point, a periodic orbit, or a heteroclinic cycle 
(or a heteroclinic network, i.e., the union of heteroclinic cycles). 

A proof of this theorem can be found in [26]. 
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